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ABSTRACT 

The "drifting" subpulses exhibited by some radio pulsars have fascinated both ob- 
servers and theorists for 30 years, and have been widely regarded as one of the most 
critical and potentially insightful aspects of their emission. Moreover, Ruderman & 
Sutherland (1975), in their classic model, suggested that such regular modulation was 
00 \ produced by a system of subbeams, rotating around the magnetic axis under the action 

of ExB drift. Such "drift" sequences have thus been thoroughly studied in a number 
' of pulsars, but it has proven difficult to verify the rotating-subbeam hypothesis, and 

thus to establish an illuminating connection between the phenomenon and the actual 
physics of the emission. 

Here, we report on detailed studies of pulsar B0943+10, whose nearly coherent 
sequences of "drifting" subpulses have permitted us to identify their origin as a system 
of subbeams that appear to circulate around the star's magnetic axis. We introduce 
several new techniques of analysis, and we find that both the primary and secondary 
features in the star's fluctuation spectra are aliases of their actual values. We have 
£^ ■ also developed a method of tracing the underlying pattern responsible for the observed 

^ \ sequences, using a "cartographic" transform and its inverse, permitting us to study 

C$ • the characteristics of the polar-cap emission "map" and to confirm that such a "map" 

in turn represents the observed sequence. We apply these techniques to the study 
of three different Arecibo observations: a 1992 430-MHz sequence which includes a 
■ transition from the star's highly organized "B" profile mode to its disorganized "Q" 

I— i ' mode; a 1972 430-MHz "B"-mode sequence; and a 1990 111-MHz "B"-mode sequence. 

The "B"-mode sequences are consistent in revealing that the emission pattern 
consists of 20 subbeams, which rotate around the magnetic axis in about 37 periods 
or 41 seconds. Even in the "Q" mode sequence, we find evidence of a compatible 
circulation time. The similarity of the subbeam patterns at different radio frequencies 
strongly suggests that the radiation is produced within a set of columns, which extend 
from close to the stellar surface up though the emission region and reflect some manner 
of a "seeding" phenomenon at their base. The subbeam emission is then tied neither to 
the stellar surface nor to the field. While the origin of the "memory" responsible for 
the stability of the pattern over several circulation times is unknown, the hollow-conal 
form of the average pattern is almost certainly the origin of the conal beam forms 
observed in most pulsars. 

Key words: MHD — plasmas — pulsars: general, individual (B0943+10) — radia- 
tion mechanism: nonthermal — polarisation 



I. INTRODUCTION 

We wish to reopen several questions which have lain dor- 
mant for most of 20 years: What is responsible for the spec- 
tacular sequences of "drifting" subpulses observed in many 
pulsars with conal single (Sa) profiles? What is the topology 



and polarisation of the individual beaming elements which 
correspond to these subpulses — and how do they accrue to 
form the average profile? How can we understand a range of 
other observed effects — profile and polarisation modes, mul- 
tiple P3 values, "absorption" and subpulse "memory" — in 
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terms of fundamental emission processes and propagation 
effects within the pulsar magnetosphere? 

From the initial identification of "drifting subpulses" by 
Drake & Craft (1968), many suspected the phenomenon to 
be an unusually clear manifestation of the physical processes 
behind pulsar emission. Fluctuation-spectral analysis (taken 
over from scintillation studies) found its place almost imme- 
diately in assessing pulse sequences for periodicity (Lovelace 
& Craft 1968), and strenuous efforts were made by several 
groups to study the subpulse modulation in the then known 
stars (i.e., Taylor et al 1969; Lang 1969; Cole 1970; Sutton et 
al 1970; Slee & Mulhall 1970). Backer's (1970a-c,1971,1973) 
thesis research carried these early efforts to understand drift- 
ing subpulses to a new level; he systematized the then ex- 
isting studies, clarified how sightline geometry could pro- 
duce distinct modulation characteristics, developed the cru- 
cial technique of applying fluctuation-spectral analysis to 
each narrow longitude interval within the emission window 
(and of following the varying phase of a feature with lon- 
gitude), and introduced the terms within which the "drift" 
phenomenon is now almost universally understood. 

Of course, the very existence of individual subpulses 
raises unresolved questions about how these distinct ele- 
ments of emission are produced — but stars with periodic 
subpulse modulation (drifting in single stars with single pro- 
files and longitude-stationary modulation in those with dou- 
ble ones) have suggested, from near the time of their discov- 
ery, a system of regularly spaced subbeams which rotate 
progressively around the magnetic axis of the star (Rud- 
erman 1972). Several efforts have been made to delineate 
the beam topology and polarisation characteristics of these 
subbeams, but none have thus far been very successful. Pul- 
sar B0809+74 has attracted special attention because its 
subpulse modulation is so remarkably regular (Sutton et al 
1970; Taylor & Huguenin 1971) — giving a strong, high-"Q" 
feature at some 0.090 cycles/period in its fluctuation spec- 
trum. However, the most creative efforts at measuring the 
polarisation of its individual-pulse sequences have resulted 
in only marginal signal strength relative to the noise level 
(S/N), so that these characteristics are now known only in 
the average of several drift sequences (Taylor et al 1971). 

Pulsar B0943+10 is another star in 0809+74's class. Its 
fluctuation spectra also exhibit narrow features — this time 
at about 0.47 cycles/period — but its unusually steep spec- 
trum makes the star difficult to observe on a single-pulse 
basis at frequencies above 300 MHz. Its pulse-to-pulse mod- 
ulation is then nearly odd-even, with alternate pulses ap- 
pearing to progress slowly from the trailing to leading edge 
of the window over 15 or so periods. No entirely consis- 
tent picture has emerged from the existing studies of its 
drift phenomenon. Taylor & Huguenin (1971) first associ- 
ated the drift pattern with another, weaker fluctuation fea- 
ture at 0.065 cycles/period and attributed the stronger fea- 
ture mentioned above to the nearly "alternate pulse mod- 
ulation". Backer (1973) noted the general possibility that 
such a primary feature was an alias of the actual fluctuation 
frequency, but found no means of distinguishing between the 
possibilities. Backer et al (1975) then showed that the two 
features have an exact harmonicity if the secondary one is 
the aliased second harmonic of the primary one, but then 
went on to discuss the phase function of the primary fea- 
ture without considering whether it also might be aliased. 



Finally, Sieber & Oster (1975) clearly discuss some of the 
aliasing possibilities, attempt to distinguish between them, 
but (in our view) come to the wrong conclusion — and thus 
find no cause to mention any harmonicity between the fea- 
tures. We note that these various studies measured signif- 
icantly different frequencies for the primary feature, which 
are summarized in Table I. 

PSR 0943+10 exhibits two distinct profile "modes" at 
102 MHz, designated "B" and "Q" by their discoverers, Su- 
leymanova & Izvekova (1984). In an earlier paper (Suley- 
manova et al 1998), the properties of these modal sequences 
are explored at 430 MHz for the first time using an ex- 
ceptionally strong, 18-minute polarimetric observation. The 
contrasting characteristics of the two modes are striking. 
The B (for "bright") mode entails a steady, highly orga- 
nized pattern of subpulse emission — the drifting subpulses 
discussed above — which is confined to the early part of the 
emission window; whereas the weaker "Q" (for "quiescent") 
mode produces disorganized, but occasionally much more 
intense, subpulses throughout the entire window. The pri- 
mary and secondary polarisation modes (PPM and SPM, 
respectively) are emitted in both B- and Q-mode sequences, 
the former well dominated by the PPM and the latter only 
slightly dominated by the SPM, respectively — resulting in 
different levels of aggregate linear depolarisation. 

These circumstances raise many questions about how 
the pulsar's average, polarised modal profiles are constituted 
of their respective pulse sequences, which were only partially 
addressed in the earlier paper. Furthermore, these subpulses 
directly reflect the electrodynamic processes responsible for 
pulsar emission — that is, (according to the received "car- 
toon") primary particle acceleration associated with local 
"sparking" regions near the polar cap, which in turn gen- 
erate a secondary plasma "bunches" that radiate coherent 
radio emission. We will attempt to relate our observations 
and analysis to the general picture elaborated by Ruderman 
& Sutherland (1975; hereafter R&S). 

Finally, something is known about 0943+10's overall 
emission geometry. Estimates of its magnetic latitude an- 
gle q and impact angle f3 were made by one of us (Rankin 
1993a, 6) on the basis of the best information then available. 
It was clear even from these initial values that our sightline 
barely grazes the pulsar's emission cone, reconfirming its 
classification as having a conal single (Sd) profile geometry. 
Moreover, at 400 MHz and above, the pulsar's decreasing 
conal beam radius is apparently responsible for the pulsar's 
extremely steep spectral decline (Cornelia 1971). Only re- 
cently have weak detections been reported at frequencies 
above 600 MHz: Deshpande et al (1999) at 840 MHz and 
Weisberg et al (1999) at 21 cms. At and below 100 MHz 
the pulsar is one of the brightest in the sky; indeed, to our 
knowledge it is the first and still the only pulsar discov- 
ered at a frequency below 300 MHz (Vitkevitch et al 1969), 
and it is unique in exhibiting no spectral turnover down to 
some 25 MHz (Izvekova et al 1981; Deshpande & Radhakr- 
ishnan 1992). Q §11 discusses our observations, and § §111— IV 
the star's modulation features and aliasing. §VI assesses the 
star's emission geometry. §VII discusses our subbeam imag- 

* Malofeev (1999) now has evidence suggesting that the pulsar's 
spectrum begins to turn over between 34 and 25 MHz. 
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ing technique and applies it to the 1992 "B"-mode sequence. 
The "Q"-mode sequence is discussed in § VIII. §IX discusses 
the polarization- mode structure of the "B"-mode sequence. 
§§X-XI discuss the older (1974) 430-MHz sequence and a 
more recent 111.5-MHz one. Finally, §XII briefly summa- 
rizes the arguments and conclusions P] 



II. OBSERVATIONS 

The single-pulse observations used in our analysis below 
come from two programs carried out at the Arecibo Ob- 
servatory over a long period of time. The newer 430-MHz 
observation was made at the Arecibo Observatory on the 
19th October 1992 and is identical to that considered in the 
companion paper (Suleymanova et al 1998). Use of 10-MHz 
bandwidth, across which 32 channels were synthesized by 
the 40-MHz Correlator, and a 1006 [is dump time reduced 
dispersion delay across the bandpass to negligible levels. The 
resolution was then essentially the dump time of 0.328°. The 
observational procedures will be described in a forthcoming 
paper (Rankin, Rathnasree & Xilouris 2000). 

The older 430-MHz observation was carried out on the 
2nd January 1972 with a single-channel polarimeter of 2- 
MHz bandwidth and 1-ms integration time, giving a nominal 
time resolution of about 1.1° longitude. The polarimetry 
scheme is described in Rankin et al (1975). 

The 111.5-MHz observation was made on the 17th Jan- 
uary 1990 in an earlier phase of the 40-MHz-Correlator- 
based polarimetry program. Use of 2.5-MHz bandwidth and 
retention of 64 lags (channels) along with a 1406 [is dump 
time resulted in an effective resolution of about 1.265°. No 
continuum source observation is available to calibrate Stokes 
parameter V for this observation, so the baseline levels, 
which were dominated by galactic background noise, were 
thus used to calibrate the relative gains of the channels. 
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Figure 1. Longitude-resolved fluctuation spectrum for pulsar 
0943+10 at 430 MHz. A 256-point fit was used and averaged over 
the first 816 pulses of the 19 October 1992 observation. The body 
of the figure gives the amplitude of the features, with contours 
at intervals of 0.012 mjy 2 , which reach a maximum value of 0.11 
mjy 2 . The average profile (Stokes parameter /) is plotted in the 
left-hand panel, and the integral spectrum is given at the bottom 
of the figure. The pulse-sequence amplitude was adjusted before 
analysis to remove the slow secular decrease [which is intrinsic 
to the pulsar, see Suleymanova et al (1998)], so that any small 
1// "tail" in the spectrum has been suppressed. Note that four 
distinct features are clearly visible in the integral spectrum, the 
two major features as well as a pair of weak features on either 
side of the principal one (see text). 



III. PRIMARY FLUCTUATION FEATURES 

We first compute nominal 256-point fluctuation spectra for 
0943+10's "B" mode, using the 816-pulse sequence at 430 
MHz discussed in Suleymanova et al, and the result is shown 
in Figure |lj The slow (about 30%) decrease of intensity over 
this interval (which they noted) was flattened before com- 
puting the spectra to avoid possible smearing of the other- 
wise high-Q features. As expected, these longitude-resolved 
fluctuation spectra (hereafter, "LRF spectra") show a strong 
feature near 0.46 cycles/rotation period (hereafter, c/Pi), 
which most fully modulates the power on the leading edge of 
the profile, but whose effect can be discerned throughout the 
"B"-mode profile. The remarkable strength and narrowness 
of this feature is evident in the integral spectrum, from which 
its frequency can be accurately determined as 0.465 ±0.001 
c/Pi. This value agrees well with that determined previ- 
ously by Backer et al (1975), but is incompatible with those 
of Taylor & Huguenin (1971) and Sieber & Oster (1975) 
(which, paradoxically, are more accurately determined and 

t Some early results of this paper appear in Deshpande & Rankin 
(1999), Deshpande (1999), and Rankin & Deshpande (1999). 



consistent with each other); see Table I. The half width of 
the integral feature is less than 0.002 c/Pi, making its "Q" 
more than 200! No known pulsar, including 0809+74, has 
such a remarkably stable modulation feature. 

Nonetheless, the principal feature in Fig. |l] is not re- 
solved with a Fourier transform of length 256. Longer FFTs 
yield a progressively narrower feature, and even for a trans- 
form of length 816, the principle feature is only partly re- 
solved. On this basis, we obtain an fs value of 0.4645+0.0003 
c/Pi, which implies a "Q" greater than 500. Subpulse mod- 
ulation with such remarkable stability is unprecedented. In- 
deed, we cannot yet say whether this behaviour is typical 
or unusual even for 0943+10, but during this 14.9-minute 
"B"-mode sequence, its subpulse modulation has a bell-like 
quality — and we see in Figure [| that the sequence is undis- 
turbed by even a single null pulse. 

A second feature is seen near 0.07 c/Pi, and it modu- 
lates the emission early in the B-mode profile just as did the 
primary feature above. We calculate the frequency of this 
smaller feature as 0.0710 ± 0.0007 c/Pi. We note also that 
this feature's width is larger, about twice that of the princi- 
pal one. An FFT of length 512 yields the slightly improved 
(mean) value of 0.0697 ± 0.0005 c/Pi — and clearly resolves 
the width of the feature. 
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Table I. Modulation Feature Frequencies 



Source 


/ P (c/Pi) 


fp alias 




harmonic? 


Taylor &Huguenin 


0.477 ± 0.003 


not cons'd 


0.065 ± 0.005 


no 


Backer et al 


0.461 ± 0.010 


not cons'd 


0.078 ± 0.020 


exact 


Sieber & Oster 


0.473 ± 0.002 


0.527 ±0.002 


0.07 ±0.01 


not cons'd 


1992/430 MHz 


0.4645 ± 0.0003 


0.5355 ± 0.0003 


0.0710 ± 0.0007 


exact 


1972/430 MHz 


0.4591 ± 0.0009 


0.5409 ± 0.0009 


not determined 




1990/111 MHz 


0.4688 ± 0.0005 


0.5312 ± 0.0005 


not determined 
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Figure 2. Pulse energy as a function of pulse number of the 
816 pulses comprising the "B"-mode sequence. Notice the slight 
secular decrease in intensity, its periodic tendencies, and, most of 
all, notice that there is not a single "null" pulse! 

The possible relationship between the two spectral fea- 
tures warrants detailed discussion. Of course, it is possi- 
ble to speculate, as did Sieber & Oster, that one or both 
of the features are aliases of fluctuations with frequencies 
greater than 0.5 cycles/period. Indeed, we can then con- 
firm that if the principal feature is the first-order alias 
(1 - ,f 3 = 0.4645 ± 0.0003 c/Pi) of a fluctuation whose 
actual frequency f$ is 0.5355 ± 0.0003 c/Pi, then its sec- 
ond harmonic f 3 = 2/ 3 would fall at 1.0710 ± 0.0006 c/Pi; 
and this fluctuation would then have a second-order alias at 
(falalias] = fi - 1) at 0.0710 ± 0.0007 c/Pi. If, however, 
the secondary feature appears as a first-order alias, the two 
might still have an harmonic relationship, because f' A would 
be 0.9290± 0.0006 c/Pi, which would then have a first-order 
alias (I-/3 =) 0.0710±0.0008 c/Pi. The measured frequen- 
cies accurately support the premise of an harmonic relation- 
ship as long as one or both are aliased as above. Or, as 
Backer et al concluded, it was "possible to avoid the con- 
clusion of an harmonic relationship only by assuming that 
the relationship is fortuitous". In fact, even higher order 
aliases in suitable combinations would be consistent with 
the premise of an harmonic relationship. Therefore, we can 
be certain about the harmonicity of the two features only if 
we establish their alias orders or — equivalently — the "true" 
frequencies of the features, a desired information in any case. 

In order to do so, we need to sample the fluctuations 
faster than just once a period. This may be possible by ap- 



pealing to the fact that the fluctuations are sampled also 
within the finite width of the pulse. In other words, the fluc- 
tuation spectra at different longitudes can be combined ap- 
propriately based on their longitude separations to "unfold" 
the contributions that are otherwise aliased within the 0-0.5 
c/Pi band. In practice, we can achieve this simply by Fourier 
transforming the entire time sequence, which can be recon- 
structed using the available pulse sequence and filling the 
unsampled (off-pulse) regions with zeros. Such (power) spec- 
tra can be computed for pulse sequences in suitable blocks 
and then averaged. Care must be taken to ensure that the 
sampled region is wide enough to include the entire pulse- 
emission window [i.e. untruncated) and that "baseline" is 
removed correctly, so that the spectral features correspond- 
ing to the pulsar signal are not distorted. This effective win- 
dowing of the time sequence, where only the off-pulse con- 
tributions are forced to zero, leads to a convolution of the 
spectrum of the raw (unwindowed) time sequence with the 
spectrum of the periodic windowing function. However, it 
is easy to show that this convolution does not modify the 
spectral contribution of the pulsar signal, smoothing only 
the noise and hence, in fact, enhancing the signal-to-noise 
ratio (hereafter, "S/N") of the spectral features of interest. 

The results of such a computation are given in Figure 
where only one half of the symmetric power spectrum is 
plotted as a function of frequency, normalized to the pulsar 
rotation frequency Pi. There is much to see in this diagram. 
First, note the amplitude of the frequency components at 
integral multiples of the pulsar's rotation frequency. These 
are the Fourier components of its "average" profile, which 
include contributions from its "base" profile — or what is un- 
changing from pulse to pulse — as well as from its fluctua- 
tions seen on an average. The harmonics decline monoton- 
ically to a 1/e amplitude by about harmonic number 15; 
note that this corresponds approximately to the longitude 
width of 0943±10's profile. Now notice the spectral compo- 
nents at about half-integral frequencies (actually iV± 0.535 
c/Pi), which grow steadily and peak around harmonic num- 
ber 35 — a scale which corresponds approximately to the 
width of the individual drifting subpulses, and more directly 
to P2, the longitude separation between adjacent subpulses. 
Finally, notice the near-integral frequency components (ac- 
tually N ± 0.071 c/Pi) which begin to be distinguishable 
from the precisely integral components at harmonic num- 
bers above 50 and which peak at about harmonic 70; these 
represent the Fourier components of the secondary feature. 
The relationship (i.e. the relative amplitude and phase) be- 
tween these sets of modulation components define the aver- 
age "shape" of the "drifting" subpulses, while that among 
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Figure 3. An unfolded fluctuation spectrum, obtained by contin- 
uously sampling the same sequence as in Fig. |jj Three successive, 
overlapping, 256-pulse Fourier transforms (each of length about 
2.8 X 10 5 points, when the unsampled intrapulse region was inter- 
polated with zeroes) were averaged together. There are three main 
contributions to this diagram: a) Harmonics at integral multiples 
of the fundamental 1/Pi pulsar-rotation frequency, which are re- 
lated to the Fourier components of the "average" profile. These 
are strongest at low harmonic numbers and have declined to a 
relative amplitude of 1/e by about number 15. b) Components 
falling at about half-integral frequencies [(TV + 0.535)/Pi], which 
peak at about harmonic 35, and c) a set of components at frequen- 
cies near (but not exactly at) harmonics of 1/Pi [(iV + 0.071)/Pi], 
which peak near number 70. The latter two sets of components 
represent the harmonics of the highly periodic fluctuations asso- 
ciated with the pulsar's subpulse modulation. 



the components within a set define the "envelope" followed 
by the subpulse as it "drifts" in longitude. 

Just to make things clearer, we present the fluctua- 
tion spectrum of Fig. rain a different manner by stacking 
successive sections of 1/Pi width one above the other as 
shown in Figure ^ We see immediately that the principal 
feature now falls at a frequency of about 0.535 c/Pi. The 
power at a frequency of 1/Pi and its harmonics — which com- 
prise the "average" profile of the pulsar — are shown in the 
right-hand panel. The fluctuation power at other frequen- 
cies, as a function of 1/ Pi -frequency interval, is shown in 
the body of the figure, and the integral of this fluctuation 
power, summed into a single 1/Pi band, is given in the bot- 
tom panel. This spectral representation, which we will refer 
to as an "harmonic-resolved fluctuation spectrum" (here- 
after "HRF spectrum"), can also distinguish between the 
two types of modulation encountered, namely amplitude- 
and phase-modulation. For example, the summed spectrum 
(as in the bottom panel) would exhibit nearly symmetrical 
modulation sidebands in the case of amplitude modulation, 
while a clear asymmetry will be apparent for phase modu- 
lation, as is presently the case.^J So, although the two mod- 

t Were, for instance, our sightlinc to make a central traverse 
across 0943+10's polar cap, rather than the highly tangen- 
tial one it does, we would expect to see the same primary 
modulation feature. However, it would appear as a longitude- 
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Figure 4. An unfolded Fourier spectrum as in Fig. |i[ The am- 
plitude of the frequency components at the fundamental rota- 
tion frequency 1/Pi and its harmonics are essentially the Fourier 
components of the "average" profile; and these are plotted in the 
left-hand panel of the diagram. The body of the figure gives the 
amplitude of all the other frequency components in the spectrum 
up to 100/Pi as a contour plot (with contours at amplitude in- 
tervals of 0.040 mjy 2 which reach a maximum value of 0.3530 
mjy 2 ), and the bottom panel shows the sum of these frequency 
components, collapsed onto a 1/Pi interval. The principal feature 
now falls at its true (unaliased; see text) frequency of about 0.535 
cycles/period with its first harmonic at about 1.071 cycles/period. 
We also see two "sidelobes" associated with the principal feature 
as well as a slight "leakthrough" at its alias frequency (from the 
unplotted negative-frequency portion of the spectrum). 



ulation types are in general mixed, the asymmetry between 
the 0.535 and 0.465 c/Pi features in Fig. ^ is so great that 
we can regard them as the signature of an almost "pure" 
phase modulation. 

Fig. ^presents exactly the same information as in Fig. ^, 
but in a slightly more direct way. Note that nothing has been 
suppressed here. The two figures give the amplitude of all 
the Fourier components of the PS, and its power falls pre- 
dominantly in one of the three sets of components discussed 
above. Remarkably, the aggregate power in the non-integral 
components is comparable to that in the integral ones. The 
Fourier components of the noise are not discernible in ei- 
ther Figs. |^ or ^, implying that its power (as expected) is 
distributed uniformly over fluctuation frequency and is thus 
completely negligible in any given frequency interval Again, 



stationary modulation — that is, as an amplitude modulation 
without "drift" — and its fluctuation power would thus be divided 
between positive- and negative-frequency features corresponding 
to "drift" to the right and to the left. Such pairs of identical 
features at positive and negative frequencies reliably identify the 
presence of amplitude modulation. 

§ Figs. |i| and ^ also exhibit the exceptional quality of the Arecibo 
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Figure 5. Fluctuation-phase and amplitude for the primary 
(0.535 c/Pi — dashed curve) and secondary (1.070 c/Pi — dotted 
curve) features in Figs. Kj-W Fluctuation power is plotted as a 
function of longitude within the average profile in the top panels, 
whereas the fluctuation phase is given in the lower panels. 



Fig. |^ clearly demonstrates the "bell like" quality of this un- 
precedented "B-mode" sequence from pulsar 0943+10. 

We can now be certain that the principal feature in 
Fig. [j] is the alias of a fluctuation whose actual frequency 
is greater than 0.5 c/Pi and that the secondary feature is 
its second harmonic. From Figs. ^ & ^, it is clear that 
these features, seen in the LRF spectra, are indeed higher- 
order aliases of the two sets of frequency components whose 
distances from zero frequency are harmonically related — 
that is, with maximum intensities at about harmonic 35 
and 70, respectively. Using a 512-point fft and interpo- 
lating the peak positions, we determine that the two re- 
sponses have frequencies of 0.5352 ± 0.0006 c/Pi and 
fs of 1 + 0.0710 ± 0.0006 c/Pi— accurately demonstrating 
their harmonic relationship. Further, we see that the sec- 
ondary feature has a well resolved width of about 0.0006 cy- 
cles/period, more than twice that of the primary periodicity. 
Finally, given the strength and narrowness of the primary 
feature's second harmonic, it is not ridiculous to explore the 
possibility of detecting its third harmonic, and a quick cal- 
culation (f-" [alias] = 3/3 — 1) confirms that it would fall at 
some 0.6065 ± 0.0012 c/Pi, where indeed, a minor feature 
can be seen in Fig. ^j. 

One other line of evidence can be brought to bear the 
aliasing of the observed features — that is, the fluctuation 
phase rate as a function of longitude associated with each 

40-MHz-Correlator-based PS observations. Were the noise not 
virtually "white" and free of periodic contributions, their effects 
would stand out sharply in these diagrams. 



of them. Their respective phase rates can be used to assess 
the harmonicity of the features, as well as to determine the 
drift-band spacing (P2) in longitude. Figure |E| gives average 
profiles showing both the amplitude and the phase of the 
power associated with the primary and secondary features. 
The maximum rate in the former is just less than 36°/° near 
the longitude origin, which is compatible with a longitude 
interval between subpulses P2 of just over 10° — say, 10.5°. 
Note that the phase rate is less steep under the "wings" 
of the profile, as might be expected if the subpulses move 
along a curved path. The phase rate of the secondary fea- 
ture is much faster and seems to reach some 70°/° near the 
peak of the profile, where the fluctuation power also has a 
maximum. We have not reproduced here a similar figure for 
the ostensible third harmonic, which exhibits a phase rate 
nearly thrice that of the primary feature. The harmonicity 
of these phase rates argues strongly that the 0.071 and 0.607 
c/Pi features are indeed, the aliases of the second and third 
harmonics of the primary feature and that we have correctly 
identified their alias orders. This conclusion is further cor- 
roborated by the fact that 1/P2(ms) is about 7V/Pi(ms), 
where N is the peak harmonic number (35) associated with 
the primary modulation feature. 

In summary, the above analysis has assisted us in 
unfolding 0943+ 10's otherwise multiply folded fluctuation 
spectrum and ruling out various combinations of P3 (along 
with their implied drift directions). The ambiguity, however, 
between the remaining two possible combinations — namely, 
P3 < 2Pi (thus implying negative drift) and P3 > 2Pi (pos- 
itive drift) — cannot be resolved through the above analysis, 
owing to the limited sampling of the modulation (which our 
sightline allows) during each rotation of the star. 



IV. SPACING OF THE ROTATING SUBBEAMS 

Given the regularity of subpulse drift in the current 0943+10 
sequence, and the implied stability of the drift bands, an 
important geometrical issue can offer useful insight. If the 
stable subpulse drift corresponds to a circular pattern of ro- 
tating subbeams — as is envisioned by the received "cartoon" 
of pulsar emission as well as by the Ruderman & Sutherland 
1975) theory — then we should be able to estimate the an- 
gular (magnetic azimuthal) separation between the adjacent 
subbeams using either a) the observed longitude interval be- 
tween subpulses (P2) or b) the associated rotation in polari- 
sation ang le (PA)P1 Of the two, we can use the PA-rotation 
approach immediately — at least approximately. 

We are then interested in the amount of PA traverse 
within an interval of longitude corresponding to the sepa- 
ration between subpulses — that is, P2. If the subpulse drift 
does correspond, in actuality, to a pattern of subpulse beams 
which pivot about the magnetic axis, then we might expect 
that the orientation of their linear polarization (which is 
ostensibly fixed to the projected magnetic field direction) 
would correspond closely to their magnetic azimuth. 

' 0943+10's subpulses exhibit no distinct PA signature as re- 
ported for B0809+74 (Taylor et al 1971) or B2303+30 (Gil 
1992) — rather, the PA varies simply with the longitude and po- 
larisation mode. 
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Figure 6. Composite of 4 partial polarisation profiles (a,b,d,e) folded at different phases of the nominal 1.87- Pi modulation cycle 
depicted in Fig. M as well as partial profiles of both the "drifting" power (c) and the "base" (f). All were computed from the entire 
816-pulse "B"-mode sequence of Fig. hj. Two and at times three subpulses are seen in the former which "drift" progressively from later 
to earlier phases about 2-3° from plot to plot. Note that the PA difference between subpulse peaks is negative and just less than 
30° — a circumstance which is just compatible with 20 subpulse beams and an inside (poleward) sightline traverse (see text). This is very 
clearly seen in the PA behaviour of the "drifting" (c) and "base" (f) profiles, where we find PA rates of -2.7 and perhaps about -4 °/°, 
respectively. 



Figure g(a,b,d,e) gives a composite of four average 
Stokes profiles corresponding to equal intervals of phase 
within the 1.87-Pi primary modulation cycle. As expected, 
we see the average characteristics of two (or sometimes 
three) subpulses at progressively varying positions in lon- 
gitude. These figures delineate the average polarisation of 
the drift sequence, and we will return to a fuller discus- 
sion of them below. However, note that the rotation of the 
polarisation angle (PA) over an interval corresponding to 
the spacing between the centres of these averaged subpulses 
(P2) — what we call \Pi — 1& about —28° when the subpulse 
pair straddles the longitude origin (or, in the "drifting" pro- 
file, Fig. []e). This immediately provides a direct means of 
estimating the number of subbeams which are possibly ar- 
ranged along a circular ring centred on the magnetic axis. 
Neglecting spherical effects and assuming that the subpulse 



beams fall closer to the magnetic than to the rotational axis, 
we find two possibilities: if |xp 2 | > ^M°), the magnetic az- 
imuth interval between the subpulse beams rj is |xp 2 | — P2 
and an inside (negative, poleward) sightline traverse is in- 
dicated, whereas if \xp 2 \ < P2, then 77 ~ \xp% \ + P% an d 
the sightline makes an outside (positive, equatorward) tra- 

verse|. Here, r\ ~ 28° - 10° ~ 18°, and thus 360718° in- 
dicates about 20 subbeams. A more precise determination 
of this number should be possible once the overall geometry 
has been determined, or, on the basis of some other consid- 
eration (as we will see in the following section). It follows 
that the system of such subbeams would make a full rota- 

II These approximations will be justified in §VI 
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Figure 7. Fluctuation spectra as in Fig. [l] for pulses 129-384. 
The "sidebands" of the principal feature are particularly strong 
in this sequence (see text). The contour interval is 102 mjy 2 up 
to a maximum of 714 mjy 2 . 

tion around the magnetic axis in an interval of just 20 times 
Pz — that is, around 40 rotational periods of the star. 



V. THE MODULATION ON THE MODULATION 

Given the unprecedented stability of the phase modulation 
in the 0943+10 PS discussed above, we can legitimately at- 
tempt a level of analysis which would heretofore have been 
inconceivable. It is clear that any underlying pattern of ro- 
tating subbeams is stable over many times the above 40- 
Pi scale. This immediately implies that their number must 
be an integer. Also, we would expect to see a periodic fea- 
ture corresponding to, in general, the rotational cycle of 
the subbeam pattern unless all the subbeams are identi- 
cal (or utterly random) in amplitude. Such a tertiary pe- 
riodicity, if any, would manifest itself as sidebands to the 
spectral features associated with the primary phase modu- 
lation. Whether such sidebands are symmetric or not would 
depend on whether the subbeams differ from each other only 
in intensity, or in interbeam spacing or in both, and their 
narrowness would depend on how stable the subbeam pat- 
tern is in time. 

Two other minor features can be discerned in Figs. [I] 
and ^ which fall symmetrically about the primary one, and 
their change of position in Fig. Elshows that all three fea- 
tures are aliased responses in FigTTfl. These features are seen 
even more clearly in the spectra corresponding to a small 
part of the sequence (pulses 129-384) as shown in Figure Q 
Careful measurements (using the fluctuation spectra over 
only a few degrees of central longitudes where the modula- 
tion sidebands are significantly stronger) fall at frequencies 
about 0.0268 ±0.00037 c/Pi higher and lower than does the 
primary one at 0.535 c/Pi. 



The separation of the minor feature pair from the pri- 
mary is indicative of a major frequency component in the 
fluctuation power of the tertiary amplitude modulation. 
Lower frequency amplitude modulation would produce sym- 
metrical features closer to the fundamental, and nothing is 
seen in the fluctuation spectra here. Consequently, that there 
is just one obvious such pair of "simple" — that is unimodal — 
features associated with the tertiary amplitude modulation 
suggests that it is itself here tone-like. Were this tertiary 
modulation more complex, we would expect to see additional 
pairs of symmetrical features at larger (but harmonic) sep- 
arations from the primary one and, only a slight deviation 
from the tone-like modulation is discernible in Fig. [?]. 

As already noted, the "drifting" character of 0943+10's 
pulse sequences as well as the asymmetrical nature of the 
primary feature in its fluctuation spectrum are strongly in- 
dicative of a phase modulation. The narrow and symmetrical 
pair of features about the primary one in Fig. M (and also 
in Figs. ^ & ^) strongly suggest the presence of a regular, 
highly periodic, amplitude modulation of the exquisitely reg- 
ular phase modulation associated with the star's "drifting" - 
subpulse sequences. This suggests that the putative sub- 
beams differ from each other primarily in their intensities, 
and that the intensity variations are a rather smooth and 
stable function of magnetic azimuth. 

The period of this tertiary modulation is then about 
1 /0.027 Pi /c or just over 37 rotation periods. We note imme- 
diately that this, in turn, is very close to 20 times the period 
of the fundamental phase modulation, which at 1/0.5355 
Pi/c or some 1.867 rotation periods, has a 20-cycle period 
of some 37.35 ± 0.03 rotation periods. This agreement is no 
accident or coincidence; any tertiary modulation not fully 
commensurate with the primary phase modulation would ei- 
ther broaden the primary feature or distort the "sidelobes" 
around it. Given the accuracy to which the frequencies of the 
secondary and the tertiary modulations can be estimated, 
we can now conclusively rule out the possibility that P3 is 
greater than 2 rotation periods, because 1/0.4645 Pi/c is 
simply not harmonically related to 37.35 Pi. 

The period of this tertiary modulation is then about 
1/0.027 Pi/c or just over 37 rotation periods {i.e. 37.3±0.5). 
Given that the length of the sequence analysed corresponds 
to many times this tertiary cycle and that the "sidelobe" 
features have the sharpness suggestive of a remarkably sta- 
ble modulation, an harmonic relationship between the two 
periods must exist. We note immediately that tertiary mod- 
ulation cycle is indeed very close to 20 times the period 
of the fundamental phase modulation, which at 1/0.5355 
Pi/c or some 1.867 rotation periods, has a 20-cycle period 
of some 37.35 ± 0.03 rotation periods. This agreement is no 
accident or coincidence; any tertiary modulation not fully 
commensurate with the primary phase modulation would ei- 
ther broaden the primary feature or distort the "sidelobes" 
around it. Given the accuracy to which the frequencies of 
the secondary and tertiary modulations can be estimated, 
we can now conclusively rule out the possibility that P3 
is greater than 2 rotation periods, because 1/0.4645 Pi/c 
is simply not harmonically related to 37.35 Pi. The only 
commensurate value of P3 (the one smaller than 2 Pi, as 
discussed earlier) implies with certainty that the drift direc- 
tion is negative. The observed harmonic relationship then 
also allows us to now adopt a refined estimate of the tertiary 
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modulation period as simply 20 times the better-determined 
value of P3. It is important to mention here that, while the 
integral relationship between the modulation frequencies has 
resolved the longstanding drift-direction ambiguity, the ex- 
act number of subbeams (20) is in no way crucial to our 
further analysis below. 

Therefore, we believe that this analysis has resolved the 
question of how and which way pulsar 0943+10's subpulses 
drift. Sieber & Oster's (1975) fig. 1 delineates some of the 
aliasing possibilities. The "A" band represents an apparent 
drift with a P3 of some 15, and it seems to entail gradual 
subpulse motion from later to earlier phases — that is, a neg- 
ative drift. Band "B" ostensibly represents a negative drift 
associated with a value of P3 just less than 2, and Band "C" 
one with a P3 just greater than 2. 

Stated differently, Sieber & Oster's diagram clearly indi- 
cates the nearly alternate-pulse character of the basic mod- 
ulation. So, it is not straightforward to distinguish whether 
the basic drift is to the right or to the left. Furthermore, as 
these authors emphasize, the LRF spectra in Fig. |l|— and all 
others like it — have the same difficulty. Any of the features 
in such a diagram can be aliases of fluctuations at frequen- 
cies greater than the 0.5 c/Pi Nyquist frequency associated 
with sampling at a 1/Pi rate. Only by appealing to the fi- 
nite width of the pulse can this ambiguity be alleviated as 
we have demonstrated above. 

The "B" track in Sieber & Oster's plot is then the fun- 
damental oscillation associated with 0943+10's drift, and 
the "A" track its first harmonic. This establishes the drift di- 
rection as negative. The fundamental 0.535 c/Pi periodicity 
implies an associated P3 value of some 1.867 periods — that 
is, in two rotations we see the next subpulse at a slightly ear- 
lier phase. Also, the effective "P3" of the second harmonic is 
1/1.07 or about 0.93 Pi/c — whose difference from 1.0 Pi/c 
describes the "A" track, while the conspicuously missing al- 
ternate subpulses are due to the fundamental modulation. 

The basic character of this drift behaviour can now be 
unambiguously represented by folding the entire 816-pulse 
"B"-mode sequence at the 1.867- Pi interval corresponding 
to the primary 0.535 c/Pi feature. We give such a display 
in Figure ^, where the subpulse intensity as a function of 
longitude and phase within the modulation cycle is shown 
as a contour plot in the central panel. The varying energy in 
the modulation is given in the right-hand panel and the non- 
fluctuating "base" (minimum-intensity) profile is shown in 
the bottom panel. While this display resolves longstanding 
questions about pulsar 0943+10's "drift", we now also see 
that its overall character is remarkably complex. 

Finally, the conclusion that 0943+10's "drift" is neg- 
ative is very important to arguments we will make in the 
following sections. The actual physical direction of the sub- 
pulse "drift" observed is negative — that is, it is apparently 
in the same direction as the pulsar's rotation. 

Having established that the tertiary modulation cycle 
is 20 times the primary P3 period of 1.867 Pi/c, the inter- 
pretation is clear. We therefore conclude that the observed 
"drifting" behaviour is a manifestation of the rotation of 
a relatively stable system of 20 subbeams of systematically 
differing intensity, spaced uniformly in magnetic azimuth. 

One means of verifying this circumstance is to fold the 
pulse sequence at the period of the tertiary modulation, 
and it would be prudent to use the value of this period 
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Figure 8. Block-averaged pulse sequence which corresponds to 
the principle "drift" pattern in pulsar 0943+10's "B" mode. The 
816 individual pulses are folded at the 1.867-rotation period inter- 
val, which in turn corresponds to the primary modulation feature 
(alias resolved) at some 0.5355 c/Pi in Fig. 2. The strength of 
the feature as a function of longitude and phase within the 1.87- 
period cycle is given as a contour plot in the body of the figure, 
and its integrated energy as a function of phase in the right-hand 
panel, whereas the aperiodically-fluctuating "base" (which has 
been removed from the other displays prior to plotting) is given 
in the bottom panel. Contours fall at about 26.7 mjy up to a 
maximum of 187 mjy. 

based on the phase modulation, which we calculated above. 
The result for the 816-pulse sequence under consideration is 
shown in Figure and it bears detailed discussion. The av- 
erage amplitude as a function of longitude over the tertiary- 
modulation cycle is depicted in the body of the figure; the 
integral over longitude is given at the left, and the ampli- 
tude of the aperiodically-fluctuating "base" (which has been 
subtracted before plotting the other panels) is shown in the 
bottom panel.^j 

First, note the remarkable circumstance that folding 
at the tertiary-modulation period reveals an orderly sub- 
sequence of 20 replicating emission elements, each of which 
have their own distinct relative amplitude and behaviour as 
a function of longitude. As is expected, one of the 20 ele- 
ments is significantly stronger than the others by upwards 
of a factor of two. Although each extends over 10° or more 

** We note here that this "base" profile (see also Fig. ^(f ) has a 
width of some 11-12°. It is not at all clear what interpretation 
should be made of this "steady" (as against the phase- modulated) 
contribution to the pulse sequence, but, were we to interpret this 
"base" feature as a core component, it is interesting both that 
its form is symmetrical about the longitude origin (whereas, the 
"drifting" power is not at all so) and its width is suggestive of a 
magnetic colatitude of some 11°-12° [Rankin (1990): eq.(5)]. 
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Figure 9. Block-averaged pulse sequence corresponding to the 
430-MHz "B"-mode observation in Figs. [M. The 816 individual 
pulses are folded at the tertiary modulation period, which was cal- 
culated to be 37.346 rotation periods. Average pulse intensity as a 
function of longitude is represented in the body of the figure as a 
contour plot (contours at 51 mjy up to 457 mjy). The longitude- 
integrated intensity over the period of the tertiary-modulation 
cycle is plotted in the left-hand panel, and the amplitude of the 
aperiodically-fluctuating "base" — which is subtracted from the 
other displays — is given in the bottom panel. 



in longitude, some peak much earlier than others and ex- 
hibit a surprising range of individual characteristics; the 
brightest one, for instance, has a double amplitude struc- 
ture. We emphasize that these elements must be relatively 
stable over this 15-minute "B"-mode sequence; it is their 
consistent brightness, and systematic variation in bright- 
ness, which apparently produces the amplitude modulation 
at a frequency 1/20-th that of the primary phase modula- 
tion. We do not usually think of the processes associated 
with subpulse "drift" as having stability on such long time 
scales — and it is very possible that they usually do not. 

Nonetheless, given the regularity of subpulse drift in the 
current 0943+10 sequence and the evidence that it entails a 
further modulation indicative of a repeating pattern of just 
20 elements, we are forced to reconsider received ideas about 
the "drifting subpulse phenomenon" and its possible origin 
in a pattern of slowly rotating subpulse beams (Ruderman & 
Sutherland 1975). Our results provide unexpectedly strong 
support for the simple "cartoon" understanding of "drifting" 
subpulses as a regular pattern of emission centres rotating 
around the pulsar's magnetic axis. 



ical beams associated with pulsar radiation. And, uniquely 
for 0943+10, our sightline apparently grazes its beam so 
narrowly that we miss ever more of its emission at higher 
frequencies. The pulsar is thus an excellent example of the 
conal single (Sa) class (Rankin 19936), and its value of 
B12/P 2 , which is about unity, supports this classification. 
Questions, nonetheless, remain about its geometry: Do we 
observe an inner or an outer cone (e.g., Rankin 1993a)? Does 
our sightline pass poleward or equatorward of the magnetic 
axis (e.g., Narayan & Vivekanand 1982)? Can we understand 
the pulsar's steep spectrum in terms of the conal narrowing 
with frequency? And finally, is the overall conal geometry 
consistent with the 20-beam "carousel" pattern indicated 
by the tertiary modulation? 

Information on the pulsar's geometry comes from con- 
sidering its PA traverse in conjunction with its profile width. 
Then, we can model these values using the relations in 
Rankin (1993a), eqs. (1-6). Reference to Suleymanova et 
al's fig. 4 — which treats the same 1992 430-MHz sequence — 
shows the PA behaviour for the 816 "B"-mode and 170 "Q"- 
mode pulses, and it is clear that the two modes exhibit some- 
what different PA sweep rates. Their Table II gives values 
of -2.4°/° and -3.6°/° for the "B'-' and "Q"-mode PA 
rates at 430 MHz and -3.4° /° and -3.6°/° at 102.5 MHz, 
respectively. Also important is their fig. 8, which shows the 
evolution of the "B"- and "Q"-mode profiles between 25 and 
430 MHz. A further set of time-aligned profiles from Arecibo 
observations is given by Hankins, Rankin & Eilek (2000). 

Our initial problem for 0943+10 is that we do not know 
where we cross its 430-MHz beam — that is, where j3 (the 
impact angle of the sightline with respect to the magnetic 
axis) falls with respect to p (the outside, half-power radius 
of the conal emission beam) . At higher frequencies our sight- 
line cuts well outside of the radial maximum and ultimately 
outside the half-power point p, whereas at lower frequen- 
cies the double profile forms indicate cuts inside the radial 
maximum point. Therefore, the observed dimensions of the 
430-MHz profile will probably not correspond to those given 
for the overall pulsar population by Rankin (1993a) — even if 
we could extrapolate properly to 1 GHz — because those are 
based on more central sight-line traverses that do cross the 
radial maximum point on the hollow conical emission beam. 

Let us then model the 0943+10 geometry for each of 
four possible configurations (inside traverse, inner cone; in- 
side traverse, outer cone; outside traverse, inner cone; and 
outside traverse, outer cone) as follows: a) Calculate the 
conal beam radii p outer (inner) throughout the entire frequency 
range using Mitra & Deshpande's (1999) expression [pat- 
terned after Thorsett's (1991) profile-width relation] with 
trial values of the index a 

(1 + 0.066,/ogJPr 1 m 

P Poutcr(inncr) ^ 066 ' 

where the 1-GHz values of the conal beam radii /9 ou tcr(mncr) 
are taken according to Rankin's eq. (5)0. b) Least -squares 



VI. PULSAR 0943+10'S EMISSION GEOMETRY 

We expect a pulsar with drifting subpulses to have an emis- 
sion geometry characterized by a highly tangential traverse 
through (or, perhaps better, along) one of the hollow con- 



tt While most pulsars exhibit what Rankin (1993a, b) defined as 
"inner" and "outer" cones, there is now evidence that a few stars 
exhibit cones with both larger and smaller characteristic radii (Gil 
et al 1993; Kramer et al 1994; Rankin & Rathnasree 1997; Mitra 
& Deshpande 1999). 
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fit the observed low frequency profile widths (111.5 MHz 
and below) to calculated values from the usual spherical 
geometric relation 



Af/> = 4 sin 



i(p/2 + /3/2)sin(p/2-/3/2) 



sin a sin(a + (3) 



(2) 



(Gil 1981; Rankin 1993a) with trial values of the magnetic 
co-latitude a and /3. c) Least-squares fit the frequency- 
dependent part of eq. (1) to obtain the index a, using it 
in step a) above in an iterative solution. With this proce- 
dure, one is in a position to compute the values of /3/p even 
at frequencies where f3 exceeds p. 

Two further useful quantities can be calculated here. 
One is the magnetic azimuth angle r\ that corresponds to a 
particular value of Pi . Using very simple spherical geometry, 
this angle is given (for small angles) by 

T) = 2sin" 1 [sin(P 2 /2)sin(a + /3)/|sin/3|], (3) 

where P2 is measured as a longitude interval (not a time), 
and, for 0943+10 at 430 MHz has a value of just over 10°. 
Reference to Fig. |^ gives P2 values in the 8.5-9.0° range, but 
a more accurate value comes from the fluctuation-spectral 
phase rate of the primary (0.535 c/Pi) feature in Fig. |^ — 
around 0° longitude the phase rate is some 34 ± 1°/°, which 
implies a P2 value of 10.6 ± 0.6°. [The smaller values im- 
plied by Fig. ^ probably reflect the curvature of the sub- 
beam track relative to the sightline. Note that (see just be- 
low) 10.6° cos x p 2 ~ 9°]- Clearly, if the drifting represents a 
pattern of 20 subbeams around the periphery of the polar 
cap, then this angle 77 must take a value of just 18° . 

Another geometrical parameter associated with the 
drifting pattern is the PA rotation between subpulses — what 
we called "PA P2" above. This PA rotation directly reflects 
the relative curvatures of the sightline and the conal beam. 
Its value is easily calculated from the usual single-vector 
model (Radhakrishnan & Cooke 1969; Komesaroff 1970) by 
evaluating X p 2 = %{ p< i/ 2 ) ~ x(-^/2), or 



X p 2 = 2 tan x [- 



sin a sin(P2/2) 



(4) 



"sin(f) cos a — cos(£) sin a cos(P2/2) 

where ( — a + f3. Incidentally, eqs.(3) and (4) reduce nicely 
to rationalize the rules we used in §IV. For small angles 
eq.(4) reduces to \p 2 ~ RpaP 2 (where the PA sweep rate 
R PA = sina/sin/3), and eq.(3) to 77 ~ {\Rpa\ + W/3)P 2 ; 
thus 77 ~ xp 2 + sgn/3P2 and will always be positive as long 
as the subbeams do not fall closer to the rotation axis than 
to the magnetic axis. 

Computations for all four emission geometry models 
were carried out to determine the best-fitting values of a, 
(3, and a, such that \ 2 ~ — w c )/Aw] 2 , where 

AW was taken as 0.4°. The first of these models — our pre- 
ferred one, for an inner traverse and inner cone — is given 
in Table II. The columns giving p, 13/ p, and the emission 
height h [from the centre of the star; see Rankin's (1993a) 
eq.(6)] depend only on the fitted values of a, f3, and a. w and 
w c are the respective measured (taken, wherever possible in 
the "B" mode for consistency) and modeled outside, half- 
power profile widths [the latter from eq.(2)]. p c (not shown) 
was computed for comparison from Rankin's eq.(4), and the 
quantities Rpa, r\, and \p 2 evaluated from the best-fitting 
parameters. 

A quick scan of Table II shows immediately that the 



inner cone, inner traverse model yields a very good fit and 
reasonable parameters. This is also true for the outer cone, 
inner traverse model which gives almost identical results, 
but with a 15.39° and (3 -5.69°. The two models have vir- 
tually identical values of (3/p, Rpa, rj, xp 2 , an d X 2 i so t ne 
fitting does not distinguish between them; indeed, it appears 
that models with any reasonable 1-GHz p value will behave 
similarly. 

By contrast, the outer traverse models are of two dif- 
ferent types: Inner and outer cone models with values of a 
and negative (3 similar to those above have nearly identical 
values of Rpa and \p 2 (and a about 0.53), but yield slightly 
poorer fits {\ 2 ~ 2.2). These, however, result in very differ- 
ent values of 77 — about 42° — as might have been expected 
from eq.(3), which depends strongly on the sign of (3, unlike 
the various other geometrical relations [i.e., eqs.(2) or (4)]. 
Finally, there are a pair of inner and outer cone "pole in 
cone" models with a about 2.9 and 4.0°, respectively, and 
Rpa 0.7, 77 18°, xp 2 about +7.4°, and \ 2 some 1.6. Clearly, 
none of these models are acceptable because the former fail 
to match the observed values of 77 and the latter those of 
Rpa and \p 2 ■ 

Table II. Inner Cone, Inner Traverse Emission Model 



f(GHz) 


PC) 


P/P 


h(km) 


w(°) 


w c {°) 


1.000 


4.13 


-1.04 


125 






0.430 


4.23 


-1.01 


131 


7? 




0.306 


4.29 


-1.00 


134 




0.1 


0.111 


4.49 


-0.96 


147 


16.5 


16.7 


0.103 


4.51 


-0.95 


149 


17.5 


17.5 


0.061 


4.65 


-0.92 


158 


23 


22.7 


0.049 


4.72 


-0.91 


163 


25 


24.8 


0.040 


4.79 


-0.89 


168 


27 


26.9 


0.034 


4.85 


-0.88 


172 


28 


28.6 


0.025 


4.98 


-0.86 


181 


32 


31.8 



Observational constants: P2, 10.5° 
Fitted parameters: a, 11.58°; (3, -4.29; a, 0.396; \ 2 , 1-12 
Derived parameters: Rpa, -2.7°/°; 77, 18.0°; xp 2 i -27.9° 

Our slightly circuitous means of computing pulsar 
0943+10's geometry — using only the profiles at 111 MHz 
and below with resolved double forms indicative of \(3\/p 
well less than unity — permits us to assess how the sightline 
cuts the pulsar's conal beam at higher frequencies. The row 
in Table II just above 111 MHz gives the frequency for which 
w c is near zero — that is, the frequency at which the model 
p would be just equal to \(3\. At 430 MHz, the table shows 
that our sightline passes just outside the half-power radius 
on the conal beam, by 0.07° and 0.09°, or just less than 
2% according to the two inner-traverse models, respectively. 
No wonder that 0943+10 is so difficult to detect above 400 
MHz! 

We can now draw some remarkable conclusions. We 
noted in §111 that 0943+10's subpulse drift is negative — 
that is, in the same direction as its rotation. Then, we con- 
cluded just above that our sightline makes an inside (nega- 
tive) traverse — that is, crossing between the magnetic and 
rotational axes, with a negative rotation of the PA. If we 
imagine that we are looking at the pulsar with the "closer" 
rotational pole "up", the emission beams sweep past our 
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sightline from right to left, thus implying that from our per- 
spective, the pulsar is rotating clockwise P|. 

Further, the subpulse beams are also rotating from right 
to left, but around the magnetic axis, which is "below" us; 
therefore, they rotate around the polar cap counterclockwise 
(in the frame of the star). This drift in 0943+10 is opposite 
to the Ruderman & Sutherland (1975) sense (see their fig. 
4); they discuss Ferraro's Theorem in an Appendix but do 
not clearly state its implication that polar cap plasma must 
rotate in the same direction as the star as viewed by an 
inertial observer. Ruderman (1976) gives a correct and clear 
discussion of this point. 



VII. IMAGING THE EMISSION ZONE 

The remarkable order implied by Fig. ^ of just 20 subpulse 
emission beams rotating about the star's magnetic axis de- 
mands a further level of analysis. We have therefore devel- 
oped a coordinate transformation between the usual system 
of pulsar (rotational) colatitude ( (= a + (3) and longitude 
if and a system rotating about the magnetic axis described 
by colatitude R and azimuth O. The emission beams rotate 
around the pole with a period P3 = NP3, where (just as in 
Ruderman & Sutherland 1975) N is the number of subbeams 
and P3 the usual subpulse-drift interval. 

If, further, we number the pulses by k from a reference 
pulse ko and measure the longitude y> relative to an origin 
defined by the longitude of the magnetic axis ipo, then 
is the sum of a rotation f? rot and a transformation # t rans as 
follows 



Base intensity 
-10 10 



Polar Emission Map 



e 



+ ±0trs 



(5) 



where the sign of (9 ro t will always be negative according to 
Ferraro's theorem, and the sign of # tr ans is positive for cw 
rotation (of the rotational pole which is closest to the sight- 
line) and negative for ccw rotation. Then 

e tot = 2-n:[k - k + (ip - ip )/2ir]/P 3 , (6) 

^ti-ans = sin _1 [sin^sin(v3 - <p )/ sinP] , (7) 

where 

R= 2sin _1 [sin 2 ({^-^o}/2)sinasinC + sin 2 (/3/2)]' . (8) 

Correct identification of the longitude of the magnetic 
axis ipo is very important to carrying out the above transfor- 
mation. Usually, the magnetic axis is expected to fall close 
to a profile's center, but for 0943+10 at 430 MHz, we have 
seen in the companion paper that the "B"-mode profile is 
truncated on its trailing edge, compared with its "Q"-mode 
counterpart. It is, however, only the "B"-mode profile that 
is truncated at 430 MHz, its "Q"-mode counterpart is much 
broader; if we take the center of the latter as close to the 
magnetic axis, note that this point falls close to the trailing 

tt Our discussion here is not fully consistent with the conventions 
of the rotating-vector model, which are reviewed nicely by Everett 
& Weisberg (2000). In our case, the "closer" rotational pole points 
away from the direction of the star's angular momentum vector, 
so that qhvm is 180° - a, and /3rvm is — 
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Figure 10. An image of the accessible emission zone recon- 
structed using the "cartographic" transform as defined by eqs.(5- 
9). This image is based on the observed "B"-mode sequence of 
816 pulses. The image or polar map is projected onto the po- 
lar cap shown in the main panel, with the "closer" rotational 
axis at the top of the diagram. The bottom and the side panels 
show the average- and the "base" -intensity profiles, respectively, 
as functions of the angular distance from the magnetic axis (i.e., 
magnetic colatitide). Contours indicate the sightline path at the 
top of the diagram. Here, the star rotates clockwise, causing the 
sightline to cut the counter-clockwise-rotating subbeam pattern 
from left to right. 



half-power point of the 430-MHz, 'B"-mode profile (see Su- 
leymanova et at figs. 4 & 8) . [The inflection point of the PA 
traverse should lie close to the magnetic axis (apart from 
small relativistic effects, Blaskiewicz et al 1991), but that 
point is difficult to determine accurately for this pulsar.] Fi- 
nally, we can appeal to the time-aligned average profiles in 
Hankins et al (2000). Note that some of the 430-MHz profiles 
are narrow and early (and thus presumably "B"-mode dom- 
inated); whereas another one is more extended and appears 
to be primarily a "Q"-mode average. Overall, we see here 
that the point on the 430-MHz, "B"-mode profiles which 
aligns best with the centers of their low frequency counter- 
parts is the half-power point on their trailing edge. We will 
then take this point, provisionally, as our best estimate of 
the longitude of the magnetic axis. 

With all the necessary ingredients at hand, we pro- 
ceed to perform the "cartographic" transform defined by the 
above equations (and the inside traverse, inner cone model in 
Table II) and to reconstruct an image of the accessible emis- 
sion zone, based on the observed "B"-mode sequence of 816 
pulses. The result is given in Figure [j^, with the emission- 
pattern image or "map" projected on to the polar cap shown 
in the main panel. The bottom and the side panels show the 
average- and the "base" -intensity profiles, respectively, as 
functions of the angular distance from the magnetic axis 
(i.e., magnetic colatitide). The polar map shows the distri- 
bution of the 20 nearly uniformly spaced subbeams, some 
noticeably more intense than others. The nearly tangential 
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character of our sightline traverse is evident from the ap- 
parent shapes of the individual features, which are almost 
certainly depicting only a small outer portion of the "true" 
subbeams structure. The sharply rising radial intensity dis- 
tribution (bottom panel) and the absence of any hint of a 
peak (or an inflection) over the accessible range of radii is 
also indicative of this circumstance. 

Note that we have adopted, without loss of generality, 
a simple but useful convention in this transform, wherein 
the rotational longitude associated with the magnetic axis 
defines a vertical line through the image centre, with its 
top pointing to the "nearer" rotational pole. So, a negative 
value of /3, as in the present case, would correspond to the 
observer's sightline sampling along a curved track through 
the upper half of the image. Then, the star's rotation is 
cw when the impact angle (3 has the same sign as the PA 
traverse (with respect to the observed longitude), and ccw 
when they do not. Thus, for the present case, a) the star's 
rotation is in the clockwise direction, and b) the subbeam 
pattern rotates a the ccw direction, given that the observed 
drift is towards the leading edge of the pulse window. 

At this point, it is appropriate to ask the obvious ques- 
tion: How unique or even accurate are the results given 
by this "cartographic" transformation? The uniqueness — 
or rather the correctness — of the polar map depends only, 
and directly, on the correctness of the emission geometry 
and the other parameters which are required to carry out 
the transform. It is very sensitive to the circulation time P3, 
the longitude of the magnetic axis cpo), and the sense of the 
sightline traverse (i.e., the sign of 0) — and if these quanti- 
ties are incorrect, the image will be smeared or distorted. 
The actual values of a and (3, however, largely just scale the 
image, so that their correct specification is much less crucial. 

In practice, some of these inputs that define the trans- 
form may not be known with the desired accuracy, making 
this transform less conclusive. However, this difficulty can 
be overcome by invoking an "inverse" transform, wherein 
the map reconstructed from the original pulse train is, in 
turn, used to produce a new sequence of single pulses. If the 
parameters used in the mapping transformation are indeed 
correct, then the artificial sequence should match the original 
one in its detailed fluctuation properties. Comparison of the 
two sequences is best effected on a longitude-to-longitude ba- 
sis, such that the new sequence at a given longitude is com- 
pared with the original one pulse-for-pulse and sample-for- 
sample. A simple cross-correlation coefficient, appropriately 
normalized, provides an adequate quantitative measure for 
such comparison — though, of course, those for each range 
of longitude must be combined through a suitably weighted 
average — indeed, quite a robust one, given all the possible 
ways by which the two sequences can differ. 

The inverse transform thus provides a powerful "clo- 
sure" path to verify and refine the "input" geometry. The 
inverse transform, of course, is defined by the same set of 
relations as for the (forward) "cartographic" transform. We 
have used this closure path to confirm our present assump- 
tions and also to conclusively rule out other geometries that 
were originally considered as plausible. We illustrate this 
with a simple example of two model geometries, differing 
only in their sign of (3. Shown in Figure [n]a,6 are cross- 
correlation maps of the fluctuations at different longitudes 
across the pulse window. Plotted in the left and bottom 
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Figure 12. Chaotic individual-pulse behaviour during the 1992 
"Q"-mode sequence (pulses 817-986). The average profile is 
shown in the bottom panel and the pulse-energy variation in the 
left-hand panel. Note the extremely intense subpulses and the en- 
hanced emission on the tailing edge of the window as compared 
with the "B" mode. 



panels are the average profiles for the new and the original 
sequences, respectively. It is easy to show that the diagonal 
elements of such correlation maps alone would suffice to as- 
sess the match. The other details in the displays do carry 
important information and can provide useful clues for im- 
proving the match. 

However, the most significant potential uses of the in- 
verse transform may come through studying what features 
of the observed sequence are lost in the process of recreat- 
ing the artificial sequence from the polar-cap map. Note, for 
instance, that any pattern of subbeams will produce a se- 
quence that will average to a symmetrical profile — e.g., see 
the left-hand panels of Fig. [tl| These differ noticeably from 
the asymmetric 430-MHz average profiles observed. 



VIII. THE "Q"-MODE PULSE SEQUENCE 

The 816-pulse sequence thus far considered corresponds to 
the pulsar's "B" mode; this is not the case, however, for the 
170 pulses that follow it. Suleymanova et al (1998) studied 
the attributes of these two modal sequences, and they find 
both a sharp boundary between "B"- and "Q"-mode charac- 
teristics at pulse 816/817 and evidence for slower variations 
that both anticipate this "mode change" and follow it. Here, 
we will first examine the fluctuation properties of this short 
"Q"-mode sequence, which is displayed in Figure 0— the 
full sequence, average, and pulse-energy variations are given 
in the main, bottom, and left panels, respectively. We see 
again that this "Q"-mode sequence is hardly "quiescent"; 
while being weaker overall, some individual subpulses are 
much brighter than any observed in the "B" mode. 
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Figure 11. Cross-correlations between the observed pulse sequence and an artificial one computed from the polar map using the inverse 
cartographic transform. The central panels display the correlation between different longitudes across the pulse window. Two cases are 
shown which differ only in the sign of the impact angle /3, illustrating the usefulness of the inverse transform "closure path" to assess and 
refine the geometrical parameters on which the subbeam imaging depends. The left-hand and bottom panels give the average profiles 
corresponding to the artificial and the observed sequences, respectively. The contour intervals in both plots is about 0.12 mjy 2 . 
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Figure 13. Harmonic-resolved spectrum for the 128-pulse "Q"- 
mode sequence as in Fig. [l2| The intensities have been smoothed 
to reduce the dominance of the ultra-strong subpulses. The spec- 
trum clearly indicates that the likely true modulation frequency is 
about 0.77 c/Pi, which would alias to 0.23 c/Pi in the longitude- 
resolved spectra. Contour levels are about 0.015 mjy 2 . 



The intensities of this 128-pulse sub-sequence have been 
smoothed to reduce the effect of the ultra-strong pulses. 
Its profile is noticeably broader, particularly on the trailing 
edge, and its LRF spectra show prominent modulation at an 



apparent frequency of about 0.23 c/P\. The HRF spectrum 
shown in Figure [13 clearly indicates that the actual modula- 
tion frequency is likely to be about 0.77 c/Pi, which would 
have a first-order alias at 0.23 c/Pi in the HRF spectra. 

Note that the harmonic amplitudes peak at very low 
numbers — about 1 or 2. If these sidebands are interpreted 
as a phase modulation, then the corresponding P2 would def- 
initely be in excess of 180° — that is, it would be larger than 
the pulse window and comparable to the rotation period. 
Therefore, it cannot be a phase modulation, so we should 
attribute this feature to amplitude modulation and inter- 
pret the asymmetry in the sidebands as result of a some- 
what faster modulation (with a frequency of 0.77 c/P\ or 
1.77 (2-0.23) c/Pi). We will assume the frequency is 0.77 
c/Pi and attempt to verify whether this is correct. Figure |l3| 
shows the sequence as folded at this interval, with the "base" 
profile (bottom panel) first subtracted from the main panel 
display. The modulation is deep and confined to a -2° to 
10° longitude interval. We have examined the two halves of 
this sequence see if the behaviour is similar, finding that the 
modulation over the wider part of the pulse is common to 
both subsets, while the feature at about —3° longitude is 
present only in the first half and is probably due to an odd 
strong pulse. This amplitude modulation feature manifests 
itself in the pulse-energy variations, where we see a tendency 
for every fourth pulse to be more intense than the average 
(see the left panel of Fig. [12]) . 

Apart from this prominent feature, a weak modulation 
is also visible near (about 5% less than) the primary "B"- 
mode feature discussed above, suggesting that the "B"-mode 
modulation is not entirely absent. A striking difference be- 
tween the "B"- and "Q"-mode behaviour is the longitude 
range over which the periodic modulation dominates. Later 
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Figure 14. Folded sequence corresponding to the 128-pulse "Q"- 
mode subsequence in Fig. [l2| The data have been averaged over 
10 bins spanning the amplitude-modulation feature at 1.293 Pi/c. 
The "base" profile (bottom panel) has been subtracted from the 
folded sequence displayed (main panel). The modulation is deep 
and confined to a -2° to +10° longitude range. Contours are at 
13 mjy intervals. 



longitudes are weakly modulated in the "B" mode, whereas 
just the opposite occurs here. Moreover, in the one case we 
see a virtually pure phase modulation, while here we have a 
strong amplitude modulation. Such marked modulation dif- 
ferences would be more understandable if, for example, the 
"Q" mode were associated with more central longitudes than 
those of the "B"-mode, and, while there is some indication 
that this could be the case, an error in the central longitude 
could produce much the same effect. Given all these differ- 
ences as well as the polarisation differences, it is surprising 
that both modes appear to have a nearly identical "base" 
profile — where the "base" is defined as that power which 
does not participate in the dominant systematic/periodic 
fluctuations. 

The abrupt change in fluctuation features at the "B"- 
to "Q"-mode transition appears to begin with pulse 817, 
which is both broad and very intense. However, the gradual 
decrease in "B"-mode intensity that precedes the onset of 
the "Q" mode (see Suleymanova et al 1998) is unlikely to 
be a coincidence. Hence, we have taken a closer look at the 
behaviour of the sequence just before and after the mode 
change to see if the ostensibly sudden transition could have 
been expected, or whether it is actually as sudden as it ap- 
pears. We therefore examined the sequence prior to the "Q"- 
mode onset to see if the 0.77 c/Pi feature was present, even 
at a low level — and as far as we can determine, it is not. 

In an effort to better characterize this modal behaviour, 
we asked whether the cartographic-transform images span- 
ning the "B"- to "Q"-mode transition might provide some 
useful insight. An immediate question to assess, however, 
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Figure 15. (a) Autocorrelation function (ACF) of the 128-pulse 
"Q"-mode subsequence shown of fig. [l^. The ACF shown here 
has been smoothened by a 5-lag window to reduce the corruga- 
tion contributed by the prominent periodicity discussed above, 
(b) Results of a similar analysis of "B"-mode pulses 129-384 for 
comparison. Note the peaks at a lag of about 40 periods in both 
the cases. 



was whether or not the parameters defining the transform 
for the "B"-mode sequence would be appropriate for the 
"Q" mode. Clearly, the geometrical parameters could not 
change, but we had no a priori means of knowing whether 
the subbeam rotation rate would remain fixed. While we 
were able to determine the round-trip circulation time for 
the "B"-mode regime, the "Q"-mode sequence exhibits no 
corresponding higher-order modulation on which to base an 
independent determination — as would be the case, for in- 
stance, if the subbeam pattern was not stable over a full 
rotation. In order to see whether the "Q"-mode pattern is 
stable over at least one circulation time, and, if so, to es- 
timate this interval, we examined the autocorrelation func- 
tion (hereafter, "ACF") of the pulse-energy sequence. The 
"Q"-mode ACF (solid curve) shown in Figure [13 has been 
smoothed by a 5-sample window to reduce the corrugation 
of the 0.77 c/Pi modulation. For comparison, the dashed 
curve gives a similar ACF for the "B"-mode pulses 129-384. 
The peaks at lags of about 35-40 periods for the "B"-mode 
sequence are expected; however, it is remarkable that they 
appear to persist in the "Q"-mode sequence as well. Indeed, 
that both ACFs exhibit correlation at some 35-40 periods 
provides very strong evidence that the circulation time of 
the subbeam pattern has remained essentially unchanged 
throughout the transition to the"Q" mode — despite the dra- 
matic change in modulation properties. 

Further, we applied the cartographic transform to a 
large number of short sequences, each with a large fractional 
overlap. Viewing the resulting images, in a slow "movie" - 
like fashion, allows us to "see" the evolution of the subbeam 
pattern with time. Since each of the images uses a rela- 
tively small number of pulses, the sampling is thus some- 
what sparse on the periphery of the polar cap. Thus, we have 
smoothed the image suitably to reduce this patchiness, en- 
suring that the basic details are not lost. In the early part of 
the "B"-mode sequence, most of the 20 subbeams are bright 
and nearly uniformly spaced. With time, the subbeams wane 
and wax in intensity on a scale of several circulation times, 
and overall they gradually decline in intensity. Some either 
bifurcate or partially merge with their neighbors, but soon 
reestablish their characteristic, regularly spaced, vigesimal 
configuration. Toward the end of the "B" -mode sequence, 
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most of subbeams are weak, with only one or two retain- 
ing their brightness. 

At the onset of the "Q" mode, most of the subbeams are 
weak but still distinguishable; however, first one and then 
two are exceptionally strong, with their intense discharges 
filling the area around them and spilling into larger radii 
or heights. Their lifetimes appear shorter than Pi, not al- 
lowing adequate sampling, thus resulting in their "streaky" 
character. These dramatic "Q"-mode events do not entirely 
displace the old ("B"-mode) subbeams, some of which crowd 
together, leaving space in the ring. In the process, more sub- 
beams appear, whose configuration is necessarily more com- 
pact, which leads to faster fluctuations. Groups of closely 
spaced subbeams are evident, with spacings which corre- 
spond to the observed 0.77 c/Pi feature. It would appear 
that large intervals of weak activity prompt the occurrence 
of the intense discharges. It may then be that eventually the 
closely spaced subbeams evolve so as to have just enough 
space for 20 — thus reestablishing the conditions for the 
bright and stable configuration that characterizes the "B"- 
mode sequences. 



IX. POLARISATION OF THE SUBBEAMS 

Both our analyses above and those of Suleymanova et al 
(1998) raise fundamental questions about subbeam polari- 
sation: a) How is the average profile comprised of primary- 
, secondary-polarisation mode (hereafter, PPM and SPM), 
and perhaps unpolarised (hereafter, UP) power, and what 
is the relation between the linear and circular polarisation 
(hereafter, LP and CP) in the former? b) How can we ac- 
count for the longitude offset between the PPM and SPM 
partial average profiles? c) Do the two modal profiles have 
distinct PA sweep rates, and how do these combine in the 
total "B"-mode average profile? d) How does the depolari- 
sation of the star's radiation occur, and can we understand 
it as the incoherent mixing of two fully (but orthogonally) 
polarised basis modes? e) How does the polarisation of the 
rotating beam system differ from that of the constant "base" 
profile? And, finally, f) what do the polarisation "maps" of 
the subbeams tell us about the conditions of their emission? 

In order to explore these questions, we have devel- 
oped several "mode separation" techniques. Heretofore, such 
methods have been used to segregate a pulse sequence, sam- 
ple by sample, into a pair of modal partial profiles (and 
sometimes a residual profile) [Cordes et al 1978; Gil et al 
1981; Rankin 1988; Rankin & Rathnasree 1995, 1997] ,0 but 
for our present purposes we require techniques for separat- 
ing the original sequence into a pair of polarised modal se- 
quences (and perhaps an unpolarised sequence). The first 
method, which we have called "modal repolarisation" as- 
sumes that all the observed linear and circular depolarisa- 
tion results from incoherent orthogonal-mode mixing. Under 
such an assumption, the linear and circular depolarisation 
occurs in strict proportion, as (apart from an overall angle) 
the linear can also be regarded as simply positive or negative 

§§ McKinnon & Stinebring (1998) have carried out a very inter- 
esting analysis, which also shows that incoherent mode-mixing 
can account for the depolarisation in B2020+28. 



(indeed, as it is represented on the Poincare sphere). The ob- 
ject is then to restore the observed (depolarised) sequence 
to a hypothetical pair of fully polarised modal sequences, 
preserving both the total power and the ellipticity of the 
original sequence. 

In practice these methods proceed by examining the to- 
tal intensity i and total linear polarisation / of the sequence, 
sample by sample, and then determining whether they re- 
spectively fall above or below a threshold toN, where t is 
a chosen value (usually taken to be 2.0) and on is the es- 
timated on-pulse rms noise level. If both the quantities are 
above the threshold, then the observed PA is checked to 
see whether it falls within ±7r/2 of the model PA. The de- 
tailed maths are given in the Appendix. We have tested this 
"modal repolarisation" method on pulsars for which its mo- 
tivating assumption is best exemplified — B 1929+ 10 being 
the canonical such pulsar (Rankin & Rathnasree 1997) — 
and it functions remarkably well both at modest S /N ratios 
and when virtually all of the total power is in one of the 
modes. Of course, we make no claim that these methods are 
precisely correct in a formal statistical sense, but they rep- 
resent a convenient and useful technique for our purposes. 

The sequences of some pulsars do not divide neatly into 
two fully polarised modal sequences, and 0943+10's provide 
one such instance. It appears that a portion of the polarised 
and unpolarized power in its sequences has a different origin, 
and so we have also found it useful to divide the original 
sequence three ways, into two fully polarised modal (PPM 
and SPM) sequences and a third unpolarised (UP) sequence. 
The operation of this "polarisation segregation" method is 
quite similar to the two-sequence "repolarisation" method 
above, and we both discuss its details in the Appendix and 
immediately apply it in our analvsis and discussion below. 

The three panels of Figure give the respective "B"- 
mode average profiles corresponding to the PPM, SPM and 
UP sequences. The virtually complete, elliptically orthogo- 
nal polarisation of the modal averages — and the negligible 
levels of polarisation in the latter average — demonstrate the 
efficacy of our segregation technique. Note that the PPM is 
about five times greater than the SPM, both in power and 
in peak intensity, as the effective widths of their respective 
profiles are about the same (6.7° vs. 7.1°). The PPM profile 
leads the SPM one by some 3°, and they are slightly left- 
(LHC) and right-circularly (RHC) polarised, respectively. 
Less than half (45%) of the total "B"-mode power is po- 
larised, but the greater UP profile width ( 9.7°) makes the 
modal and UP contributions to the total intensity roughly 
equal at most longitudes. 

The PA traverses of the PPM and SPM in Fig. |l| are, 
overall, quite linear, but the curves carry little information 
in the wings of the profiles, where the linear power is low and 
thus ever more influenced by the choice of the model PA tra- 
verse. Suleymanova et al found a very significant difference 
in the PA rates of the "B" and "Q"-mode profiles, and, given 
that the former is PPM and the latter SPM dominated, 
there is a case to be made. We find it hard to understand, 
though, how the PPM-SPM difference could be anywhere as 
large as the 50% (-2.4°/° vs. -3.6°/°) "B" mode-"Q" mode 
difference that they find. We have examined PA-frequency 
maps and the modal partial profiles above and conclude that 
the PPM might be as flat as -2.4°/° and the SPM as steep 
as -3.0° j° , but these seem to be outside limits. Even these 
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Figure 16. Polarised partial-average profiles corresponding to the PPM, SPM and UP partial sequences, respectively. Note that the 
PPM and SPM profiles are both highly polarised and elliptically orthogonal. By contrast, the UP profile exhibits a negligible level of 
polarisation. Note also that the PPM and UP sequence averages are comparable, both in power and in mean (scale) intensity, whereas the 
SPM average is about 5 times weaker. The scales are in mjy. 



small differences in modal PA rate, however, combined with 
the longitude offset between their power centres, is proba- 
bly sufficient to explain the peculiar form of the overall PA 
traverse, which is shallow at negative longitudes and then 
steepens significantly at positive ones. Note that the rate dif- 
ference implies that the angle between the two modes will 
be significantly less than 90° at +10° longitude, just where 
the overall traverse is steepest. 

Turning now to the polarisation of the subpulse "drift" , 
the PPM, SPM and UP sequences provide ideal instruments 
for identifying its characteristics. First, each of these se- 
quences can be folded at the basic 1.87- Pi phase-modulation 
cycle (P3) as we did in Fig. ^ above. Figure [l?] then 
depicts the behaviour of the PPM- and SPM-associated 
modulation — the former is given as a set of contours and 
the latter as a grey-scale. Their relationship is quite striking. 
We know the slopes of the "drift" tracks quite accurately; 
clearly, on average, the subpulses move through -10.5° lon- 
gitude (Pa) in the 1.867-Pi cycle, or some 5.62°/Pi. 

We see in Fig. [T?] that the PPM and SPM tracks are 
parallel as expected, though we could not have anticipated 
that the latter would fall almost exactly in between the for- 
mer, so that PPM and SPM subpulse "tracks" will always 
be found at a longitude interval of P2/2 or some 5.3°. The 
overall behaviour is a little more complex, however, because 
there is a significant phase difference between the PPM and 
SPM maxima within the phase modulation cycle. The SPM 
maximum not only lags the PPM one by about 0.25 Pi, but 
the respective intensity variations over the cycle appear to 
be almost reflections of each other — that is, the PPM ex- 
hibits a long slow increase and then an abrupt fall, whereas 
the SPM cycle appears to begin with a steep increase which 
then tails off gradually. 

In order to examine this behaviour in more detail, Fig- 
ure hq gives the result of folding each of the three partial 



sequences at the 1.867-Pi phase-modulation interval. Here 
we give only the sequence averages (full curves) and the 
relative amplitude of the "base" profiles inside them. This 
"base" has often been subtracted in displaying the results 
of our analyses and represents that power which is either 
constant from pulse to pulse or fluctuating at frequencies 
other than the phase-modulation rate. One can readily see 
here that most of the power associated with the PPM and 
SPM is fluctuating with the basic 1.867- Pi cycle — and there- 
fore associated both with the "drifting" subpulses and the 
rotating-subbeam system. Much of the UP power is also fluc- 
tuating and can be seen to largely follow the PPM "track" 
as expected; though nearly half of this power is not fluctu- 
ating in this way, and thus must have a different origin. It is 
the sum of these three contributions to the "unfluctuating" 
power that we encountered above as the triangular "base" 
profile in Fig. ^| — and there questioned whether this power, 
on account of its about 11° width, could be associated with 
the weak tail of a core component. 

We have also carried out a similar analysis using the two 
partial sequences generated by the "repolarisation" tech- 
nique, and the result is telling. Of course, any depolarisation 
within a sequence can be interpreted as orthogonal mode 
mixing and the characteristics of the original pair of fully 
polarised sequences then estimated. For 0943+10, however, 
the large fraction of the total power that we saw above in the 
UP sequence is divided equally between the modal sequences 
in the two-way split and then (because its fluctuations are 
correlated with those of the PPM) dominates the weaker 
SPM fluctuating power in its sequence. These circumstances 
both make this latter technique less useful in 0943+10's case 
and again suggest that some part of the UP power originates 
from a source other than the mixing of two fully polarised 
basis modes. This line of interpretation is consistent with 
the UP being, mostly, a weak signature of core emission, 
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Figure 18. Partial-average profiles for the PPM, SPM and UP partial sequences showing the relative contributions of the aperiodically- 
fiuctuating "base" within each partial sequence. Note that the "base" comprises a small part of both the PPM and SPM power, but 
constitutes nearly half of the UP power. Note also that it is quite symmetrical about 0° longitude. The "base" in each case reflects both 
that power which is steady from pulse to pulse as well as any which is fluctuating at a different rate than the phase-modulation rate. 
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Figure 17. Block-averaged pulse sequences computed from the 
respective PPM and SPM partial sequences as in Fig. |i| — the for- 
mer shown as a set of contours and the latter as a grey-scale plot. 
Note that both PPM and SPM "drift" along parallel "tracks", 
but at somewhat different phases within the overall 1.87-pcriod 
phase modulation cycle. 

but certainly does not prove that this is the case. We have 
also computed a Stokes profile corresponding to the "base" 
emission, and it will come as no surprise that this profile 
shows only slight linear polarisation with a PA following the 
SPM traverse and almost no circular polarisation. 

Of course, each of these sequences can also be mapped 
onto the polar cap using the cartographic transform, and 



the relation between the PPM and SPM beams is shown in 
Figure [H^. Here, the PPM emission configuration is shown 
as a set of contours and that of the SPM as a grey-scale 
map. The respective intensity scales are, of course, not equal; 
indeed, the SPM scale has been enhanced by a factor of 
about 5 for clarity. 

Note, first of all, that the SPM emission also forms a set 
of subbeams along nearly the same circle as that of the PPM. 
Indeed, we know that all these subbeams rotate ccw about 
the magnetic axis at the centre of the diagram and that the 
star's overall cw motion about its rotation axis (at the top 
of the diagram) causes the sightline path to traverse from 
right to left through the beams, just poleward of the mag- 
netic axis. The different longitude centres and modulation 
phase of the SPM and PPM power together determine the 
relative (azimuth) orientation of their respective subbeam 
systems. On average, the SPM subbeams lag their preced- 
ing PPM neighbors by an interval which is a little less than 
half the mean subbeam spacing. It is just this lag which is 
reflected in the delay of the SPM profile power centre rela- 
tive to that of the PPM in Fig. Indeed, we see there that 
the SPM power centre follows the PPM one by 3-4° out of 
the 10.5° spacing between subpulse "drift" bands {P2)- Note 
also that both the SPM and PPM beams vary markedly in 
mean intensity, and there is some correlation between the 
strength of an SPM beam and that of the PPM beams ad- 
jacent to it. Again, on average, the two beam systems both 
seem truncated at the absolute limit imposed by the sight- 
line traverse — and to about the same extent. 

Finally, we can locate the circular polarisation within 
the "drift" bands and within the subbeam system in a sim- 
ilar fashion. Figure ^ again shows the result of folding the 
sequence at the 1.87- Pi modulation-cycle interval, but here 
we see the total power (Stokes /) plotted as a contour map 
and the circular polarisation (Stokes V) given as a colour 
scale. It is immediately clear that positive (LHC) circular 
polarisation is associated with the total power — which in 
turn is dominated by the PPM emission; whereas bright neg- 
ative (RHC) circular polarisation is seen along the parallel 
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Figure 19. Images of the accessible emission zone as in Fig. [Io| 
The PPM is shown as a set of contours and the SPM as levels of 
grey. The scales of the two modes are not equal, the weaker SPM 
having been enhanced so as to be approximately as prominent 
as the PPM. Note that the SPM power falls between the PPM 
subbeams. Recalling that the subbeam pattern rotates ccw, and 
that the star's rotation sweeps the sightline through the top of the 
subbeam system from left to right, we see that the SPM emission 
generally lags the PPM beams by just less than half the interval 
between them (see text). 



"track" corresponding to the SPM. Note that the fractional 
circular can reach high levels at certain points in the di- 
agram for both modes, but that its average level remains 
low. Of course, the implied beam configuration of this cir- 
cular polarisation can also be viewed using the cartographic 
transformation, and Figure ^l| gives this result. Again, the 
total power is shown as a set of contours and the circular 
as a colour scale. Here also, we see LHC polarisation asso- 
ciated with the PPM beams as well as the region "outside" 
them — that is, further from the magnetic axis or perhaps at 
greater height along the same fields lines. Similarly, the RHC 
is found just between the PPM beams near the positions of 
the SPM beams in Fig. |l^ and also "outside" them. 

Overall, the results of our analysis in this section are 
quite clear. The rotating subbeam system that we encoun- 
tered in earlier sections is highly polarised in an elliptically 
orthogonal manner: the prominent 20-fold subbeam system 
is dominated by PPM emission, which also exhibits sys- 
tematic LHC polarisation. Then, interleaved with this PPM 
beam system is a similar, weaker SPM beam system which 
we find to have significant RHC polarisation. These circum- 
stances appear to suggest a set of very specific conditions 
within the emission region, and we will come back to a dis- 
cussion of their implications later. 
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Figure 20. Block-averaged pulse sequences computed from 
Stokes parameters I and V , respectively, as in Fig. ^| — the for- 
mer shown as a set of contours and the latter as a colour plot. 
Note that the PPM-dominated / "drifts" along a "track" that 
has significant LHC polarisation, while RHC lies along the adja- 
cent SPM "track". The circular polarisation can be quite strong 
at certain positions along these "tracks", but overall remains at 
a low level. 
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Figure 21. Images of the accessible emission zone as in Fig. [Io| 
Stokes / is shown as a set of contours and Stokes V as a colour 
plot. Note that the PPM-dominated / subbeams are LHC po- 
larised, whereas the SPM-associated region between them is RHC 
polarised. Each of these circularly polarised regions has a ten- 
dency to extend to "outside" regions, which may be indicative of 
higher altitudes in the emission region. 
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Figure 22. Harmonic-resolved fluctuation spectrum as in Fig. [j] 
for the 1972 430-MHz subsequence. Note that the primary mod- 
ulation feature is at 0.54 c/Pi — related as before to the aliased 
0.46-c/Pi feature in the LRF (not shown). Contours are at inter- 
vals of 28 mjy 2 up to a maximum of 198 mjy 2 . 



X. THE 1972 JANUARY 430-MHZ SEQUENCE 

In this section, we look at a sequence which was recorded at 
Arecibo in early 1972 (Backer et al 1975) in order to explore 
whether the conditions encountered in the 1992 observations 
are typical or unusual. We have assessed the quality carefully 
and have selected a subset comprised of pulses 201-968 (us- 
ing only the right-circularly polarised channel as the other 
was corrupted by 60-Hz interference) from this old 1000- 
pulse polarimetric observation. Here, both the single-pulse 
S/N and resolution are poor relative to the newer sequence — 
so we used a 3-sample smoothing to improve the S/N (at 
the cost of poorer resolution) but retained the original sam- 
pling. Here, too, we find the gradual, monotonic decline in 
the pulse energy with time, much like what was seen in the 
newer sequence. The average profile, however, has a more 
symmetric appearance, which is intrinsic and not a result of 
the applied smoothing. 

In order to examine the fluctuation properties of the se- 
quence, we computed both the LRF and HRF spectra. The 
latter is given as Figure |2^ and its dominant spectral fea- 
ture at 0.5409 ± 0.0009 c/Pi appears in the LRF at 0.46 
c/Pi (as a first-order alias as discussed in §111). These fea- 
tures correspond to the phase modulation associated with 
the rotating subbeams. The second-harmonic feature is not 
visible in these spectra, but is present when the smoothing 
is not performed. The harmonic amplitudes of the fluctua- 
tion features in the HRS are centred at harmonic numbers 
32 and around 65, implying a typical subpulse separation of 
about 11±0.5°. The sequence was also folded at the primary 
modulation period in the manner of Fig. ^| and exhibited a 
nearly identical behaviour. The drift band was clearly visi- 
ble, once the "unchanging" base contribution was removed 



from that folded over the modulation cycle. The modula- 
tion period (P3) used here corresponds to a refined value 
((1 — 1/2. 153) -1 Pi) at which the modulation pattern shows 
maximum contrast in pulse-energy variation across the mod- 
ulation cycle. Here there was also a slight suggestion that 
the subbeams cross our sightline twice during each complete 
circuit in magnetic azimuth; if so, it would imply that their 
angular distance from the magnetic axis is larger than |/3|, 
which is consistent with many subaverage profiles showing 
a "conal-double" appearance. Also, since the two crossings 
through our sightline are generally expected to be symmet- 
ric with respect to the "true" longitude origin, they provide 
a useful and independent estimate of the central longitude 
value which we use. The "base" profile in this subsequence, 
though, had a larger fractional intensity and a half-power 
width of about 10°, consistent with the "base" width in the 
newer observations, periods is compatible with a feature at 

We also computed the phase variation of the modula- 
tion feature (at fa) as a function of longitude, and the results 
were virtually identical to Backer et aPs fig. 4 (though we 
have taken care to specify the central longitude correctly). 
The rate of change of the modulation phase with respect to 
longitude gives an independent estimate of P2, which we find 
encouragingly consistent with that suggested by the HRF 
spectrum as well as with that found for the 1992 sequence. 
Based on this similarity, we estimate [using eq. (3)] the sec- 
ondary modulation period to be about 20 times P3. 

Unfortunately, the relatively noisy character of this se- 
quence has frustrated our efforts to use the intensity auto- 
correlation function to confirm our estimate for the round- 
trip circulation time. The modulated intensity here is a 
smaller fraction of the average intensity than in the 1992 
observation. Also, there is noticeably more emission at longi- 
tudes corresponding to the profile wings (i.e. at outer radii) 
where the intensity fluctuations are high, but it does not 
appear to have any definite periodicity. This "outer" emis- 
sion is clearly seen in the average profile "wings" as a pair 
of weak, roughly symmetrical features. 

We have attempted to map the average subbeam distri- 
bution of this sequence, using a circulation time of 37.014 Pi 
and the "inside, inner" geometry discussed earlier. The cen- 
tral longitude estimate was verified using a bootstrap pro- 
cedure based on the inverse transform "closure" . Figure poj 
shows the result after smoothing the map with a circular 
function having an area 8 pixels to reduce its noisy appear- 
ance. The overall character of the distribution is quite simi- 
lar to that of the 1992 "B"-mode sequence; however this map 
shows significant emission at outer radii — probably from an 
outer ring of more or less randomly spaced subbeams. This 
"outer" emission is different in kind from that of the "Q" 
mode at a similar radius, basically because its distribution 
is so "conal" but random in azimuth. 

To summarize, this "old" sequence seems to show all the 
primary "B"-mode characteristics, together with some weak 
and random emission possibly associated with the same 
sources as the "Q"-mode emission. The secular decrease in 
pulse energy may well be an indication of an intense, but 
sporadic "Q"-mode sequence to follow. The modulation fre- 
quency is a little higher than that of the recent observation, 
but it is well within the variations seen in short sections of 
the more recent sequence. 
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Figure 23. Polar emission map obtained via the cartographic 
transform for the 1972 430-MHz observation in Fig. [l^. The map 
is smoothed with a circular 8-pixel-area window to reduce the 
noise in its peripheral regions. Again, only the RHC channel was 
used in the analysis. 



XI. THE 1990 JANUARY 111.5-MHZ SEQUENCE 

Now we examine at a sequence observed at a longer ra- 
dio wavelength. It consists in more than 1000 pulses, and 
we have looked at the first 1024 pulses in detail. At this 
lower frequency the profile has two well resolved compo- 
nents separated by about 10° longitude, implying that here 
the radiation-cone radius p is larger than the impact angle 
f3. This is not at all surprising if we are observing emission 
from the same subbeams whose outer periphery we imaged 
at 430 MHz. We would expect to see them more fully at 
lower frequency, because of the ever larger emission-cone 
radius. Indeed, our favored geometrical model in Table II 
indicates a conal beam radius p some 0.26° larger than at 
430 MHz, or, said differently, f3/p is here 0.96 as opposed 
to 1.02 at the higher frequency. This is a major advantage; 
and we will then hope to "see" a larger portion of the actual 
radial extent of the subbeams. 

Let us first examine the fluctuation properties of this 
sequence. The dominant feature, even at this frequency, is 
associated with the nearly alternate-pulse modulation that 
we meet yet again. Over the 1024-pulse sequence, the fluctu- 
ation frequency varies by perhaps 4% of its mean value, but 
appears stable over intervals of a few hundred pulses. We 
show, as a typical example, detailed spectra for pulse num- 
bers 1-350. Other portions of the observation have similar 
spectra, except for small shifts in the modulation frequency. 
Figure uA shows the LRF spectra, where the aliased fea- 
ture appears at a frequency of about 0.46 c/Pi. The HRF 
spectrum (not shown) confirms that the 0.46-c/Pi feature 
is again the first-order alias of a 0.53-c/Pi modulation and 
that it is phase modulated. Its second harmonic was found 
close to 0.07 c/Pi. The harmonics related to both the rota- 
tional and the phase modulation show a distinctly different 
envelope compared to their equivalents at 430 MHz. The 
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Figure 24. Longitude-resolved spectra (using 256-point fits) for 
a subsequence (pulses 1-350) of the 1990 111.5-MHz observation 
(main panel) and the total spectrum (bottom panel) as well as 
the average profile (left panel). Note the similarity between these 
spectra and those of the 430-MHz observations. The intensity 
units are arbitrary. 



difference is a direct reflection of the double-resolved profile 
as against the single form at 430 MHz. The LRF spectrum 
confirms that the harmonics of the modulation feature peak 
at about 32 c/Pi, again implying a P2 between 10-11°. 

Next, we look at the result of folding the 350- pulse se- 
quence at the corresponding modulation period (see Fig- 
ure ^) . A prominent drift band associated with the leading 
component is seen in the left half of the main panel. Appear- 
ing as a faint continuation is the "track" associated with the 
"drift" of the weaker trailing component. The "base" profile 
subtracted from the data in the main panel is given in the 
bottom panel. Note that it here resembles the average pro- 
file, unlike what we observe at 430 MHz. It is more likely, 
therefore, that this "base" reflects aperiodically-modulated 
residual emission, stemming possibly from variations be- 
tween the drift-bands of different subbeams, rather than be- 
ing the result of another form of radiation. The variation 
in the modulation phase (not shown) is remarkably similar 
to that seen at 430 MHz, implying about the same value 
of P2 as determined earlier. The conclusions based on the 
first set of 350 pulses are seen to be valid throughout the en- 
tire sequence, apart from small differences in the modulation 
frequency. It may be worth recalling that such differences, 
though on a reduced scale, were also evident in the recent 
430-MHz sequence. 

We now proceed to perform the cartographic transfor- 
mation. Having noted that the modulation properties are 
nearly identical to those for the 430-MHz "B"-mode se- 
quence, the conclusion that the round-trip delay is quite 
close to 20 times the mean P3 is unavoidable. However, since 
the differences in the P3 values for the different sections of 
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Figure 25. Subsequence folded over the modulation cycle for the 
same 1990 111.5-MHz observation as in Fig. [24] (main panel); the 
"base" profile (bottom panel) has been removed. Note the two 
distinct tracks associated with the two components seen in the 
111.5-MHz average profile. This should be contrasted with the 
single "drift track" seen at 430 MHz (see Fig. ^). 
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Figure 26. Polar map obtained via the cartographic transform 
of the same 1990 111.5-MHz subsequence as in Fig. ^ (main 
panel). The map has been smoothed with a circular (8- pixel-area) 
window to reduce the noise in the peripheral, poorly sampled 
areas. Note the less elongated appearance of the subbeam features 
as compared to those in the 430-MHz maps, implying that the 
more interior 111-MHz sightlinc samples the beams more fully. 
The intensity scales are arbitrary. 



this data are found to be different by more than the avail- 
able resolution in P3 (which is simply the inverse of the 
number of pulse-periods in the sequence), we should per- 
form the transform on suitable sections of the data sepa- 
rately, unless we wish to view only the average distribution 
of the subbeams. With the geometry already defined, we use 
the different data sections and the inverse transform "clo- 
sure" to confirm our estimate of the central longitude, and 
then examine the subbeam distributions (in simple projec- 
tion onto the polar cap) for the different sections separately. 
Figure [26] shows one such map, again corresponding to the 
first 350 pulses. The other sections of the data show a simi- 
lar distribution of subbeams, differing only in their detailed 
intensities and the locations of individual subbeams. As an- 
ticipated, we now see a much larger fraction of each of the 
subbeams than was possible at 430 MHz. The radial profile 
(in the bottom panel) shows an intensity maximum well out- 
side the /3-cutoff. In fact, the fraction of the subbeam that 
is not visible now is quite comparable to that visible at 430 
MHz. By comparing the radial profiles at the two frequen- 
cies, we can now estimate the factor by which the radiation- 
cone size has changed — 1.03 between 430 and 111.5 MHz, 
which is in good agreement with the modelling in Table II. 



XII. SUMMARY AND CONCLUSIONS 

In this first paper of a series, we develop techniques and 
present detailed studies of the drifting subpulses of pul- 
sar B0943+10, based on three sets of Arecibo observations. 
Pulse sequences from this star are remarkable for their ex- 
traordinary stability and apparent complete lack of null 
pulses. Future papers in this series will feature sequences 
of several different kinds from other low frequency observa- 
tories as well as overall analysis and physical interpretation 
of their significance. 

Our main results may be summarized as follows: 

• We have resolved the question of whether 0943+10's 
primary spectral feature is aliased and whether the sec- 
ondary feature is its second harmonic: The primary feature 
at about 0.46 c/Pi is the first-order alias of the fundamental 
phase-modulation, which has a true frequency of about 0.53 
c/P\. The secondary feature is then indeed a second har- 
monic (at about 1.07 c/Pi), which is seen as a second-order 
alias at about 0.07 c/Pi. 

• Three new techniques have been developed and applied 
to resolve the above questions: a) A "harmonic resolved fluc- 
tuation" spectrum uses the information within the finite 
width of the pulse to achieve a Nyquist frequency of 1 c/ Pi, 
showing clearly that the primary feature is aliased. Then, 
the harmonicity of both b) the harmonic amplitudes of the 
two features and and c) their phase rates are used to argue 
that the former appears as a first-order alias and the latter 
as a second-order one. 

• On this basis, we were able to determine that P3 is some 
1/(0.53 c/Pi) or 1.87 Pi and to demonstrate this circum- 
stance by folding the entire sequence at this interval. This 
implies that that drift is negative, or, in the same direction 
as the star's rotation. 

• Using several different techniques, we were able to mea- 
sure the longitude interval between adjacent subpulses P2, 
as well as the polarisation-angle rotation corresponding to 
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this interval XP2J an d to relate these to the azimuthal angle 
between adjacent between subbeams around the magnetic 
axis of the star rj. The most reliable values of Pi come from 
either the principal harmonic number of the phase modula- 
tion or the phase rate associated with its primary feature, 
both of which yield values of about 10.5±0.5° . Interestingly, 
we find that |xp 2 |, which just less than 30°, should exceed 
P2 for an inner sightline traverse, and vice versa; and that 
77 is approximately the difference between these quantities. 
Thus r\ must be near 18° for 0943+10, implying a system of 
some 20 subbeams. 

• A pair of "sideband" features are observed in the overall 
fluctuation spectra (and very strongly in certain sections). 
We identified these features as an amplitude modulation on 
the primary phase modulation. Such features can only occur 
if this tertiary amplitude modulation is stable and harmon- 
ically commensurate with the primary phase modulation. 
The separation between the "sidebands" and the primary 
feature A/ is just 1/20 of the 0.54 c/Pi frequency, indicat- 
ing that the tertiary modulation is produced by a repeating 
pattern of just 20 elements. The overall modulation period 
is then 20 times P3 or just over 37 Pi, which was dramat- 
ically demonstrated by folding the entire sequence at this 
interval. It is at this point that we have conclusive evidence 
that the aliasing question is resolved, because 0.5355/A/ 
is an integer within its errors, whereas other possible alias 
frequencies of the primary phase modulation are not. 

• Pulsar 0943+10's drifting-subpulse pattern can be com- 
pletely understood as resulting from a system of 20 sub- 
beams rotating around its magnetic axis in an interval of 37 
rotation periods or about 41 seconds, where some of the sub- 
beams are stronger than others and maintain this intensity 
difference for several circulation times. 

• We have developed a new technique, involving a "car- 
tographic" transform and its inverse, to map and study the 
underlying subbeam pattern. The forward transform merely 
expresses the observed pulse sequence in terms of a pulsar- 
frame magnetic colatitude and azimuth, rotating at the sub- 
beam circulation rate, rather than the usual parameters of 
pulse number and phase (or longitude). Techniques are de- 
scribed whereby the forward transform can be used to study 
and display the sequence, whereas the inverse transform can 
be applied in a "search" mode to determine the geometry or 
other poorly known parameters of the sequence. 

• All the techniques including the "cartographic" trans- 
forms were applied to three different Arecibo observations, 
two at 430 from 1992 and 1972 and one at 111 MHz from 
1990, and entirely compatible results were obtained for each. 

• The polarisation characteristics of the most recent ob- 
servations were also carefully investigated. Using new meth- 
ods of segregating the observed sequences into two or 
three subsequences with definite polarisation states, we have 
mapped the pulsar's emission in the two polarisation modes 
and plotted their characteristics. The PPM is about five 
times stronger than the SPM, and both closely follow a R&C 
behaviour at every longitude across the profil e]^] — though 
they appear to have slightly different PA sweep rates. Pos- 

11 We see no evidence of the anomalous PA rotation within the 
drifting subpulses as has been reported for 0809+74 by Taylor et 
al (1971) and for 2303+30 by Gil (1992). 



itive circular polarisation is associated with the PPM and 
negative with the SPM. A part of the unpolarised power, 
which drifts, seems to be the result of modal depolarisation, 
but another part seems to have an independent origin. 

• Polarised maps of the beam configuration show clearly 
that the PPM emission consists of a 20-fold system of bright, 
well confined subbeams, whereas the SPM emission is ob- 
served as a series of "beams" interleaved between the PPM 
system as well as at larger radii (or heights). This and other 
features of the beam configuration suggest that the pul- 
sar radiates at various altitudes within a system of plasma 
columns which have their "feet" near the stellar surface. 

• Our qualitative result that the pulsar's drifting sub- 
pulse emission is produced by a rotating subbeam system 
as well as our quantitative values for both the circulation 
time and physical subbeam spacing (see also Deshpande & 
Rankin 1999 a) appear fully compatible with the Ruderman 
& Sutherland (1975) model, though our analysis is com- 
pletely independent of this model. 

• The methods of the present analysis should be useful 
for studying the origin of pulse modulation in other pulsars. 
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APPENDIX. POLARISATION-MODAL 
ANALYSIS TECHNIQUES 

Here we give the mathematical details of the polarisation- 
mode separation techniques which we used in §IX. As these 
techniques proceed sample by sample, we are always con- 
cerned with the polarisation state of each individual sample, 
which can be represented by Stokes vector s = (i, q, u, v) in 
relation to the standard deviation of all sources of noise 
in any one Stoke s parameter a. The total polarised power 
is then p = q 2 + u 2 + v 2 , and the linear power / = 

y/q 2 + u 2 , where (f> = | tan -1 (u/q) is the sample polari- 
sation angle (PA). Both techniques require a model PA tra- 
verse, which is usually taken to describe the primary polar- 
isation mode (PPM) and which, for the sample in question, 
has a value (j> m . 

Our modal "repolarisation" technique follows that out- 
lined for profile polarisation in Rankin & Rathnasree (1997); 
however, we here apply it to the individual samples of a 
pulse sequence. To this end we require a few more defi- 
nitions: The total fractional polarisation / = p/i, so, for 
the primary mode I' = l/f and v' = v/f. Then, for the 
secondary mode l" = -(I' - l)/2 and v" = ~(v" - v)/2. 
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In that the two modes are fully polarised by assumption, 
%' = \Jl> 2 + v' 2 and i" = \/l" 2 + v" 2 . Finally, q' = V cos 20 
and u' = I' sin 2(j>, just as q" = I" cos 2<p and u" = I" sin 2(f). 

Insignificant portions of the resulting sequences must be 
noise-like, just as are the natural ones. Therefore, ft here rep- 
resents a noise source — that is a Gaussian-distributed ran- 
dom variable with zero mean and a standard deviation equal 
to a/2. In generating the two repolarised modal sequences, 
four different computations are carried out depending upon 
whether the sample falls above or below two different thresh- 
olds, one pertaining to its total intensity relative to the noise 
level a and a second to determine whether it is significantly 
linearly polarised. Only when a sample has significant power 
but negligible polarisation is the model angle used to re- 
construct the modes; then we compute the necessary linear 
I* = y/i 2 — v 2 and so q* — I* cos 2</> m and u* — 1* sin2<j)m. 
The overall procedure is then as follows: 



I > 2a and p < i 
(partial polarisation) 

s P (s) = ip',q',u',v') 

Ss( P ) = (p",q",u",v") 
For \4>-4> m \ < f (> f) 

I < 2a and i > 2a 
(depolarisation) 

s p = {i/2,q"/2,u*/2,v/2) 
s s = {i/2,-q*/2,-u*/2,v/2) 



I > 2a and p > i 
(hyperpolarisation) 

s P (s) = (i,q,u,v) 
s s( P ) = (h,n,n,n) 



For \<f> - 



<?(>?) 



I < 2a and i < 2a 
(noise) 
s p = (i/2,q/2,u/2,v/2) 
s s = {i/2,q/2,u/2,v/2) 



The three-sequence polarisation segregation method is 
slightly simpler to carry out. The same thresholds are used, 
h is defined as a/3, and the overall procedure is as follows: 



I > 2a and p < i 
(partial polarisation) 

s P (s) = (p,q,u,v) 
s s ( P ) = (h,h,h,n) 
s u = (i — p, n, n, h) 
For \4>-<j> m \ < f (> f) 

I < 2a and i > 2a 
(depolarisation) 

s p , s s = (n, n, h, h) 
s u = (i,q,u,v) 



I > 2a and p > i 
(hyperpolarisation) 

s P (s) = (i,q,u,v) 
s s( P ) = (n,n,h,h) 
s u = (n, n, n, h) 
For \<t>-<t>m\ < f (> f) 

/ < 2a and i < 2a 
(noise) 

s p ,s s = (i/3,q/3,u/3,v/3) 
s u = (n, h, n, h) 
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